AD-A248  307 


On-line  point  positioning 
with  single  frame 
camera  data 


Final  report 


Th.  Kersten,  K.  R.  Holm,  A.  Gnien 
Institute  of  Geodesy  and  l^otogrammetiy 
ETH-HdnggCTbeig,  CH  -  8093  Zurich 
Switzerland 


92  4  06  169 


92-08908 


UNCLASSIFIED 

SECURITY  CLASSIUCATION  Of  THIS  PAGt 


REPORT  DOCUMENTATION  PAGE 


fotm  Approved 
0MB  No  0704  01S8 
ftp  0*tt  tun  30.  lies 


I*.  REPORT  SECURITY  CLASSIEICATION 

UNCLASSIFIED 


lb.  RESTRICTIVE  MARKINGS 


2a  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  OECLASSiflCATiON/OOWNGRAOiNG  SCHEDULE 


3  DISTRIBUTION /AVAILABILITV  OF  REPORT 

APPROVED  FOR  PUBLIC  RELEASE: 
DISTRIBUTION  UNLIMITED 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S} 


S.  MONITORING  ORGANIZATION  REPORT  NUM6ER($) 

R&D  5366-EN-Ol 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

Institute  of  Geodesy  & 
Photogrammetry 


6b.  OFFICE  SYMBOL 

(If  tpplktble) 


7a.  NAME  Of  MONITORING  ORGANIZATION 

USARDSG-UK 


6c.  ADDRESS  (Gty,  State,  tnd  HBCodt) 

ETH-Hoenggerberg,  CH-8093  Zurich 
Switzerland 


7b.  ADDRESS  (City,  State.  »nd  ZIP  Codt) 

PSC  802  BOX  15 
FPO  AE  09499-1500 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

Army  Research  Office 


Sb.  OFFICE  SYMBOL 
(If  app/icab/e) 

ARO 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

DAJA45-  89-C-0004 


Be.  ADDRESS  (City,  Stitt,  ind  ZIP  Codt) 

PO  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

ELEMENT  NO. 

NO 

NO. 

ACCESSION  NO 

311C2A 

IL161102BH 

J  01 

AR 

y\.  TJTLE  (lnc/ud«  Security  CUssificetion) 

ON-LINE  POINT  POSITIONING  WITH  SINGLE  FRAME  CAMERA  DATA 


IZ  PERSONAL  AUTHOR(S) 

PROFESSOR  A.  GRUEN 


I3a.  type  of  REPORT 

FINAL 


13b  TIME  COVERED 
FROM  (H  Cl.  ri  TO 


14  DATE  OF  REPORT  (Ytir,  Month,  Diy) 

-03  ~  'fS' 


15  PAGE  COUNT 


16.  SUPPLEMENTARY  NOTATION 


17. 

COSATI  CODES  | 

18.  SUBJECT  TERMS  (Contmut  on  rtverst  if  ntetsury  ind  idtntify  by  Woe*  number) 

FIELD 

GROUP 

SUB-GROUP 

DvfCimC  rinihaL  Vdiiit  ltiv.fii.^l<\Hoi^^ 

DSU5 

Seyfithhiit  Ohci  S'*iut/ahrc^is  (it/^u-shnneuif^  iluindtz  oftfct-h  dn. 

19.  ABSTRACT  (Continue  on  reverse  if  necessery  and  identify  by  block  number) 


20  DISTRIBUTION /availability  OF  ABSTRACT 
Q  UNCLASSlFIEOtUNLIMITED  GO  SAME  AS  RPT 


13  OTIC  USERS 


Z1  abstract  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

JERRY  C.  COMATI _ 


22b  TELEPHONE  (Includt  ArtiCodt) 

071-402-7331/8490 


22c  OFFICE  SYMBOL 

AMXSN-llK-RE 


DO  FORM  1473,84  MAR 


B3  APR  edition  majr  bf  used  until  eahausted 
AM  other  editi^t  are  obsolete 


StCURITY  CLASSlFlCATlOW  OF  THIS  PAG6 

UNCLASSIFIED 


On-line  point  positioning  with 
single  frame  camera  data 


Final  technical  report 
by 

Th.  Kersten,  K.  R.  Holm,  A.  Gruen 
March.  1992 


United  States  Army 

EUROPEAN  RESEARCH  OFHCE  OF  THE  U.S.  ARMY 
London  England 

CONTRACT  NUMBER . . 

Institute  of  Geodesy  and  niotogiammetry 
ETH-HOnggerbeig,  CH  •  8093  Ziirich 
Switzerland 


Approved  for  Public  Release;  distribution  unlimited 


Table  of  Contents 


CHAPTER  1 _  Project  data  i 


CHAPTER  2  Introduction  2 


CHAPTER  3 _  Goals  of  the  project  3 


CHAPTER  4 _  Hardware  aspects  3 


CHAPTER  5  _  Adopted  algorithms  4 

5.1  Computation  of  initial  values  4 

5.1.1  Space  resection  4 
S.U  Relative  crieatatian  6 
S.IJ  Spttiil  intenection  7 

S-1-4  Ttoee-dimcnsioiia]  similari^  mnsfonnatian  8 

S.IJ  Strategy  of  providing  initulvahiet  for  the  on-line  data  procesting  8 

Adjustment  9 

5.11  Simultaneowadjuatmentliyaioledvfactorizatianndbacksiiliatilulian  9 
311  Sequential  adjustment  by  Givens  transfonnations  11 
Quality  control  13 

31.1  Bhnder  detection  and  locadon  13 
311  Systematic  emrcompenaaliim  IS 


^APTER  6 


Software  modules  of  OLTRIS  is 


Table  of  Contents 


7.1  The  user  inter&ce  18 

7.2  Inidalizatioa  18 

7.3  Measmement  and  data  processing  18 

lA  Quality  control  21 

7.5  Documentation  21 

7.6  Results  21 

CHAPTER  8 _  Conclusions  22  _ 

CHAPTER  9 _  Acknowledgment  22 

CHAPTER  10 _  References  23 _ 

10.1  Project  references  23 

10.2  General  references  24 


On-line  pomt  posiiianiiig  with  single  frame  camera  data 


On>line  point  positioning  with  single  frame  camera  data 

Final  report 

Th.  Kenten,  K.  R.  Holm,  A.  Graen 
Institute  of  Geodesy  and  Photogrammetry 
ETH-Hdnggerbeig,  CH  -  8093  Zurich 
Switzerland 


Summary 

The  project  “On-line  point  positioning  with  single  firame  camera  data”  was 
concluded  on  December  31,  in  1991  after  a  period  of  three  years.  This  final  re¬ 
port  presents  the  investigations  and  results  of  this  project  The  goal  of  the 
project  was  to  devdop  a  data  processing  concept  for  on-line  high  accuracy 
point  positioning  system.  An  operational  on-line  system,  installed  in  a  UNIX 
environment  ur.der  the  programming  language  “C’,  for  aerial  phototriangu¬ 
lation  working  on  an  analytical  plotter  or  on  a  digital  workstation  should  pro¬ 
vide  high  speed  of  operation,  accuracy  of  results,  and  case  of  data 
management  Establish^  tests  showed  that  the  Sun  Mkrosystems  worksta¬ 
tions  have  sufficient  memory  capacity  and  speed  performance  to  run  on-line 
triangulation  applkathms.  The  adopted  algorithms  in  the  on-line  triangula¬ 
tion  system  for  the  computation  of  fcdtial  values,  for  the  sequential  and  simul- 
tanous  adjustment  and  for  quality  control  will  be  prcsoited  in  this  report  A 
description  of  the  developed  software  modules  and  the  program  will  also  be 
given.  Although  this  system  was  designed  and  developed  specifically  for  aerial 
triangulation,  it  is  anticipated  that  fhture  applications  for  on-line  triangula¬ 
tion  algorithms  and  methods  will  be  found  in  robotics  and  industrial  quality 
control. 

1.  Project  data 

The  project  has  been  defined  as  "On-line  point  positioning  with  single  frame  camera  data”. 
The  research  rqiotted  in  this  final  rqxxt  has  bwn  made  possible  throu^  die  support  and 
qxmsonhip  of  the  U.S.  Government,  Department  of  Defense,  throu^  its  European  Re¬ 
search  Office  of  the  U.S.  Army  in  London.  The  period  for  research  and  investigations  was 
fixed  at  3  years,  finxn  January  1, 1989  to  December  31, 1991.  Researchers  on  the  project 
have  been  Dr.  Knut  R.  Htrim  as  Postdoctoral  Researcher  from  April  1, 1989  to  March  31, 
1990  and  Thomas  Kersten  as  Research  Associate  ftom  May  1, 1989  to  Deconber  31, 1991. 
The  investigations,  implementatioos  and  results  of  diis  project  have  been  described  six  In¬ 
terim  Reports,  which  are  listed  in  section  10.1.  (see  Project  References).  The  fiiud  r^iort 
summarizes  all  work  undertaken  in  this  project 
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2.  Introduction 

Aerooiangulation  is  an  efficient  phtnogrammetric  method  of  point  detennination  with  laige 
H«t«  sets  and  plays  an  important  role  in  surveying,  lb  improve  the  potential  efficiency,  ac¬ 
curacy  and  relial^ty  of  off-line  aerooiangulation,  the  operator  shoi^  receive  a  continuous 
information  feedback  on  the  progress  and  quality  of  his  measurementsAvoik.  In  tins  con¬ 
text,  measurements  in  large  photogrammenic  blocks  require  quality  control  at  an  early 
stage  of  data  acquisition.  W^e,  in  general,  the  procedures  of  die  conventional  off-liiw 
aeiDtiiangulation  are  separated  in  data  preparation,  data  acquisititm,  bundle  adjustment, 
blunder  detection  and  location,  and  finally  quality  control,  the  on-line  solution  integrates  all 
these  data  processing  phases  into  one  process,  ^ause  of  the  advantages  of  the  sequential 
solution,  on-line  data  processing  has  been  an  imptHtant  issue  in  photogrammetry  for  nuiny 
years.  In  general  terms,  “on-line  triangulation”  is  the  procedure  of  measuring  immedi¬ 
ately  processing  data  for  point  positioning  purposes.  “On-line”  means  direct  communica¬ 
tion  with  a  computer  during  the  triangulation  procedure,  so  that  measurements  be 
immediately  processed  after  being  acquired. 

After  the  foundation  of  a  woddng  group  “Analytical  On-line  ’Triangulation”  (WG III/4)  at 
the  Xmth  ISP  Congress  in  Helsinki  1976  the  term  On-line  Triangulation  (OLT)  has  been 
defined  in  the  resolution  T  Ulfi  of  the  XIVth  Congress  of  ISPRS  in  Hamburg  1980,  be¬ 
cause: 

"  ...  the  on-line  capability  to  measure  and  immediately  process  data  increases 
the  speed  and  reliability  of  photogrammetric  triangulation  and  may  significant¬ 
ly  improve  the  organization  of  routinely  performed  work.” 

An  increase  of  qieed  and  reliability  of  the  overall  triangulaticHi  procedure  can  primarily  be 
achieved  by  controlling  blunders  removing  false,  or  adding  new  obsovations  at  an  ear¬ 
ly  stage.  In  order  to  perform  these  functirais,  an  OLT  system  requites  effective  hardware 
and  software  for  processing  large  data  sets  in  near  real  time.  Essentially  then,  the  perfcn- 
inance  of  such  systems  depends  on  the  efficient  tuning  of  algorithms,  hardware  arid  soft¬ 
ware.  For  this  reason,  full  utilization  of  the  potential  provided  by  the  hardware  always 
depends  on  efficient  software.  The  most  time  consuming  aspect  in  OLT  is  the  adjustment  of 
the  data.  Consequently,  fast  sequential  algorithms  have  been  suggested  for  adjusting  photo- 
gramrtKtric  measurements.  Furthermore,  die  ctmqiutation  of  initial  values  as  discu^ed  in 
section  S.  1  and  the  updating  of  the  covariance  nuitrix  of  the  residuals  (Qw)  must  be  solved 
efficiendy  (section  5.3). 

Simply  stated,  the  realizing  of  an  OLT  system  requires  powerful  cranputer  hardware  and 
optimal  algorithms. 
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3.  Goabof  the  project 

The  purpose  of  this  project  was  to  develop  a  data  processing  concept  and  a  software  inqile- 
mentaticm  for  an  on-line  high  accuracy  point  positioning  system.  To  this  end  an  operational 
on-line  system  for  aerial  phototriangulation,  whose  performance  would  considerably  ex¬ 
ceed  the  currently  available  systems  with  respect  to  speed  of  operation,  accuracy  of  results, 
and  ease  of  data  management  was  to  be  built  The  system  should  be  capable  of  monitoring 
on  an  analytical  plotter  or  on  a  digital  woilcstation  using  digital  images.  It  was  recognised 
that  investigations  into 

•  Initial  value  computations, 

•  Sequential  and  simultaneous  bundle  adjustment, 

•  Blunder  detection  and 

•  Self  calibration 

would  be  needed.  Furthermore,  the  option  to  expand  the  on-line  system  by  the  use  of  kine¬ 
matic  GPS-data  and  automatic  point  transfer  was  to  be  considered.  Later  the  system  inay  be 
modified  to  be  effectively  used  in  terrestrial  applications  (close-range  photogrammetry)  as 
well.  This  expansion  to  a  more  general  on-line  triangulation  system  requires,  in  particular, 
further  investigations  into  initial  value  computations.  These  goals  are  further  eloboiated 
upon  the  project  in  the  initial  proposal.  It  should  be  noted  that  it  was  necessary  to  change 
some  of  these  goals  or  their  pricnities  and  to  introduce  new  goals  due  to  experiences  made 
in  developing  the  on-line  triangulation  system  OLTRIS.  For  example,  ato  discussions 
with  the  project  sponsor,  a  higher  priority  was  given  to  building  an  opetaticmal  user  inter¬ 
face  based  on  a  Vfindow  System  befme  implementing  self  calibration. 

4.  Hardware  aspects 

Before  developing  and  irrqilementing  an  OLT  system,  conqruting  hardware  needs  to  be 
considered  in  order  that  acceptable  perfomaance  for  the  OLT  ^rplications  can  be  achieved. 
The  computing  environment  at  the  Institute  of  Geodesy  and  Photogrammetry  consists  of  a 
netwmic  of  Sun  Microsystems  workstations.  Consequently,  OLTRIS  has  bm  developed 
and  implemented  in  a  UNIX  environment  (SunOS  Release  4.1.1)  on  a  Sun  wmkstation  in 
the  programming  language  “C”.  It  was  assumed  that  memory  capacity  of  the  system  be 
practically  unlimited  with  respect  to  the  number  of  image  points,  object  points,  photo¬ 
graphs  and  cameras  in  the  so^are.  Moreover,  response  times  should  be  well  within  the 
limits  of  the  qrerator’s  patience.  Therefore,  at  the  conmtencement  of  this  project  the  neces¬ 
sary  efiSciency  and  ciqucity  of  the  SUN  workstation  was  evaluated.  The  possible  mraiory 
requirements  ate  given  in  the  following  exanqile.  For  a  normal  equation  system  based  cm 
otoervation  data  of  SO  photos  and  500  object  points,  each  measun^  in  4  images,  2  Mb  ate 
requited,  while  100  photos  aixl  KXX)  object  points  need  S  Mb.  This  data  resides  in  memory 
during  sequential  adjustment  Also,  for  digital  images  (or  windows  thereoO  available 
on-screen  viewing,  10  image  patches  ofS12x512at  8-Ut  resolution  pixels  occupy  2.S  Mb. 
However,  only  a  few  need  to  be  in  menxiry  at  any  giv»  time.  Furthermore,  for  digirireH 
pmnts,  which  can  be  defined  by  small  patdiN  Qess  than  25  x  25  pixels),  KXX)  image  patch¬ 
es  requite  about  0. 5  Mb.  Again,  not  all  would  be  required  to  be  in  monoty  at  any  time. 

Using  the  workstations  of  SUN  Microfystenis  diese  tequitements  of  memoy  capacity  and 
speed  att/canbe  met  The  efficiency  tff  three  generatioas  of  standard  SUN  wortetations  is 
given  in  taUe  1.  In  this  context  dw  qreed  of  die  computer  is  judged  by  "MiSion  Floating 
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Point  Operations  Per  Second"  (MFLOPS)  and  by  "Million  Instructions  Per  Second" 
(MIPS),  using  double-precision  LinPack  as  reference. 


workstation 

speed 

frequency 

capacity 

SUN  3/110 

IMFLOFS]  (MIPS) 

0.4  2.0 

IMHzl 

16.7 

(Mb]  Mwwmnf  (liA]  mgdnMJ 

12  12 

SUN  SPARCstation  1+ 

1.7 

15.8 

25.0 

16 

16 

SUN  SPARCstation  2 

4.2 

28.5 

40.0 

16 

96 

Table  1 :  Performance  of  SUN  workstations 

Specific  examples  of  the  performance  of  a  SUN  workstation  for  OLT  triplications  are  pre¬ 
sented  below  in  section  5.2.2  and  5.3. 1. 

5.  Adopted  algorithms 

Of  particular  research  interest  ate  the  algorithms  used  for  computation  and  providing  initial 
values  for  the  unknown  parameters,  for  the  procedures  of  the  adjustment  and  for  the  quality 
control  of  the  measurements.  The  algorifiuns  used  in  OLTRIS  are  described  below. 

5.1.  Computation  of  initial  values 

For  simultaneous  and  sequential  adjusunent  OLTRIS  has  to  provide  initial  values  for  the 
unknown  parameters,  e.g.  object  point  coordinates,  exterior  orientation  parameters.  The 
adopted  approaches  for  computing  initial  values  ate  summarized  in  the  following. 

5.1.1.  Space  resection 

The  spatial  position  of  the  camera  and 
the  rotation  angles  of  its  optical  axis 
during  die  moment  of  exposure  can  be 
determinated  by  space  resection.  Hius, 
space  resection  is  a  method  of  comput¬ 
ing  the  six  exterior  raientation  ele¬ 
ments  (Xo,  Yg.  9,  m,  K)  of  one 
photognqih.  Mathenutical  formual- 
tion  of  space  resection  is  based  on  the 
collinearity  equations.  The  geometric 
principle  of  these  equations  is  that  an 
objea  point  Pj,  its  image  point  P'j  and 
the  perspective  center  O  ate  coUinear 
(Figure  1). 


Figure  1:  Geometry  c/f  space  resection 
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The  collineaiity  equations  are  expressed  by  the  following  formulas 


Fy=  (y-yp)  ^ 

where 


u 

X-Xo 

V 

=  D 

Y-Yo 

w 

Z-Zo 

and: 


X,  y . image  cooidinates, 

Xp,  yp,  c . elements  of  interior  orientation. 

X,  Y,  Z . object  coordinates  of  point  Pj, 

D . 3  by  3  orthogonal  rotation  matrix  and 


Xo.Yo,Zo....  ..  object  cooidinates  of  perspective  center  O. 

The  elements  of  interior  orientation  (Xp,  yp,  c)  define  the  location  of  the  center  of  perspec¬ 
tive  in  the  image  coordinate  system  (pnysically  the  finont  nodal  point)  and  they  are  given  in 
the  camera  cidibration  data  file  of  OLTRIS.  The  exterior  orientation  elements  define  the  lo¬ 
cation  of  the  perspective  center  in  object  space  (Xp,  Yq,  Z^)  and  the  orientation  of  the  im¬ 
age  cootdinate  system  with  respect  to  the  object  coordinate  system  (three  rotations  tp,  co,  k). 
For  the  computation  of  the  exterior  orientation  parameters  at  least  three  known  object 
points  must  appear  in  the  image.  The  image  point  coordinate  observations  of  the  object 
points  yield  six  collineaiity  equations  to  compute  the  extericn  orientation  elements  without 
redundancy.  To  exploit  any  redundancy  in  the  observations,  least  squares  adjustment  is 
commonly  employed. 

Space  resection  will  not  yield  a  solution  when  all  image  points  lie  on  a  straight  line  or  too 
tightly  grouped.  To  avoid  these  cases,  the  system  tests  Ae  configuration  of  the  available 
points  in  image  space  before  the  adjustment.  The  test  criterion  used  is  the  area  formed  by 
the  object  points  in  image  space.  The  area  formed  by  three  points  will  be  computed  as  a  tri¬ 
angle  (see  Figure  2).  If  there  are  more  than  three  points  available,  die  area  will  be  computed 


Figure  2:  Area  of  a  triangle  in  image  qiace  (left)  and  of  a 
group  of  points  in  image  space  (right) 
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as  the  boundary  of  the  extreme  points  (see  Figure  2).  The  test  criterion  is  5%  of  the  area  of 
an  image,  e.g.  for  an  image  23  cm  x  23  cm  it  is  26.45  cm^.  If  the  computed  area  is  under 
this  threshold  the  resection  will  not  be  computed  and  the  system  will  suggest  to  the  opera¬ 
tor  to  measure  more  known  points  or  to  perform  a  relative  orientation  of  the  images . 

For  the  least  squares  adjustment  of  the  space  resection  starting  values  (approximations)  of 
the  unknown  parameters  are  necessary.  As  a  start,  Xq  and  are  the  X-  ^  Y-components 
of  the  gravity  center  of  the  object  coordinate  system.  The  height  of  die  projection  center 
can  be  approximated  with  the  formula 


Zo  =  Z„.  +  c^ 

O  OS  d. 


(2) 


with: 

Zg . height  of  the  projection  center, 

^ . height  of  gravity  center  in  object  coordinate  system, 

c . camera  constant, 

do . averaged  distance  in  object  space  and 

di . averaged  distance  in  iinage  space. 


Initially,  the  rotation  angles  are  set  to  q>  =  co  =  0,  with  k  computed  as 

K  =  arc  tan— -arc  tan  ^  (3) 

yu  *os 

with; 

K . rotation  angle, 

Xjg,  yj„ . coordinates  of  gravity  center  in  image  space  and 

Xog,Yos . coordinates  of  gravity  cento’  in  object  coordinate  system. 


Using  these  approximations,  the  exterior  orientation  elements  will  be  computed  in  an  least 
squares  adjustment  of  the  resection. 

5.1.2.  Relative  orientatioii 

If  the  exterior  orientation  elements  of  two  photographs  can  not  be  determinated  by  resec¬ 
tion  in  ^lace,  the  relative  position  and  attitude  of  the  images  in  a  stereoscopic  pair  with  re¬ 
spect  to  each  other  can  be  conqinted.  The  two  photographs  are  cnientated  relatively  to  each 
other  if  all  corresponding  pair  of  rays  from  the  two  images  to  the  object  points  intersect  in 
space.  Because  of  model  itefonnations,  the  exact  intersection  in  space  cannot  be  made  in  all 
corresponding  pairs  of  rays.  For  relative  orientation,  a  minimum  of  five  corresponding 
pairs  ^  rays  are  needed,  while  a  sixth  pair  is  useful  for  checking.  In  OLTRIS  two  methods, 
the  swing-swing  method  and  the  one-projector  method,  arc  used  for  relative  orientation. 
Both  metiiods  are  based  on  the  copUuutrity  equation.  Geotnetrically  the  coplanarity  equa¬ 
tion  represoits  a  plane  which  is  formed  by  two  conjugate  image  points,  the  perqiective  cen¬ 
ters  and  the  objM  point  all  form  (condition  b  in  Hgure  3).  Tb  realize  the  coplanarity 
condition,  the  volume  of  the  tetrahedra,  as  dqiicted  in  condition  a  of  figure  4,  must  be  zero. 
The  formulas  of  die  ct^lanarity  equation  are  described  in  Karan  (ed.)*. 

^  Nop-Topmi'Hihicri  Pboit^nBsniBBy»  AmviuD  SodBQf  of  Flmo^HuniMiy,  19t9,|ii|«  43/46. 
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The  swing-swing  method  determinates  two  rotation  angles  (<p,  k)  of  the  left  image  and 
three  rotation  angles  (tp,  to,  x)  of  the  right  image,  with  the  other  orientation  parameters  be¬ 
ing  fixed.  This  method  will  be  only  used  in  OLTRIS  for  confuting  the  rotation  angles  of 
the  first  model.  The  one-projector  method  determinates  the  y-  and  z-base  components,  and 
the  rotation  angles  of  the  second  image,  with  the  first  seven  orientation  parameters,  six  of 
the  first  image  and  one  of  the  second  image,  being  fixed.  If  redundant  image  point  observa¬ 
tions  are  measured,  the  orientation  is  computed  by  a  least  squares  adjustment. 


5.1.3.  Spatial  intersection 

If  the  exterior  orientation  elements  of 
two  photographs  in  a  overlapping  pair 
(stereopair)  are  known,  the  spatial  posi¬ 
tion  (X,  Y,  Z)  of  any  point  can  be  deter¬ 
minated  fiom  the  measured  photo 
coordinates  in  the  two  images  using  the 
method  of  spatial  intersection.  The  geo¬ 
metric  principle  (Ingure  4)  is  based  on  a 
spatial  intersection  of  two  rays  from  two 
fixed  points  (camoa  stations  0|,  O2). 


The  foimulas  for  spatial  intmection  used  in  OLTRIS,  are  essoitially  the  same  as  shown  for 
this  purpose  in  Karara  (ed)*. 

*  HuBiook  of  Woii-T«niogiipl»icil  Photagramiieliy,  Ammcm  SocieCy  of  Phottgimiiiwliy.  1979.pige8S-89 
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The  method  gives  initially  four  equations  for  solving  three  parameters  (X,  Y,  Z).  Those  four 
equations  can  be  combined  in  three  different  ways,  each  giving  three  linear  equations.  In 
Karaia,  one  such  combination  is  chosen,  giving  one  set  of  equations.  Then  the  three  equa¬ 
tions  are  combined,  giving  a  “closed  solution”  for  X,  Y,  and  Z,  i.e.  one  formula  for  each  pa¬ 
rameter  (X  =  ffobservations),...).  The  formulation  fails  in  the  event  of  an  “unfortunate” 
base  direction,  due  to  a  denominator  approtrehing  zero.  The  algorithm  in  OLTRIS  therfore 
contains  all  the  three  sets  of  closed  solution  formulas.  The  used  set  of  the  appropriate  for¬ 
mulas,  chosen  to  give  best  precision  in  the  result,  has  the  largest  denominator. 

The  computed  object  point  coordinates  are  located  in  the  coordinate  system  fixed  by  the  ex¬ 
terior  orientation  parameters  of  the  images. 

5.1.4.  Three-dimensional  similarity  transformation 

Because  of  insufficient  datum  information  during  the  triangulation  processing,  the  object 
points  and  the  perspective  centers  must  be  defined  in  a  local  coordinate  system.  Later,  a 
three  dimensional  coordinate  transformation  can  be  used  to  transfer  the  local  into  a  global 
datum  as  long  as  a  minimum  of  three  identical  points  in  both  coordinate  systems  are 
present,  l^thin  OLTRIS,  a  three-dimensional  similarity  transformation  transforms  a  local 
into  a  global  coordinate  system,  which  is  defined  by  the  ground  control  points.  If  insuffi¬ 
cient  control  point  information  is  available,  a  local  coordinate  system  is  fixed  at  the  start  of 
the  triangulation.  The  formulas  used  for  the  transformation  are  ffie  same  as  described  in  Al- 
bertz/Kreiling^.  The  parameters  of  the  transfOTmation  to  be  determinated  are  a  uniform 
scale  factor,  three  rotations  and  three  translations.  These  unknowns  can  be  obtained  by  a 
standard  least  squares  adjustment.  Note  that  a  minimum  of  three  points  should  be  known  in 
each  coordinate  system  -  two  X-,  Y-coordinates  and  three  Z-coordinates.  The  three-dimen¬ 
sional  similarity  transformation  provides  OLTRIS  with  the  flexibility  to  start  the  triangula¬ 
tion  in  a  local  coordinate  system  and  to  transform  it  to  the  global  system  whenever 
sufficient  control  point  information  is  available. 

5.1.5.  Strategy  of  providing  initial  values  for  on-line  data  processing 

A  special  problem  at  the  start  of  the  photogranunetric  triangulation  is  to  provide  initial  val¬ 
ues  for  the  unknown  parameters  (i.e.  object  point  coordinates  and  exterior  orientation  ele- 
.nents)  of  least  squares  adjustment  Herein,  the  choice  of  functions  for  computing  initial 
values  depends  on  the  given  information.  Therefore,  the  triangulation  procedure  is  based 
on  one  of  the  following  cases: 

•  exterior  orientation  elements  given,  e.g.  determined  by  General  Positioning  System 
(GPS),  no  control  points  (rare  case). 

•  some  control  points  given,  deteimined  by  geodetic  network  adjustment;  no  known  ex¬ 
terior  orientation  elements  given  (normal  case). 

•  control  points  and  known  exterior  orientation  elements  given  (ideal  case). 

•  no  information  from  a  ground  coordinate  system;  global  transftnmation  impossible 
(bad  case). 


^  Phatognminetric  Ouide,  Wicfaminn  Veriig,  Kwlnuhe  1980,  ffp  45  ■  47. 
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In  each  case,  a  strategy  for  providing  the  necessary  initial  values  is  needed.  For  this  reason, 
besides  the  coding  of  functions  for  die  computation  of  initial  values,  the  investigations  of 
the  project  have  also  concentrated  on  finding  the  best  strategy  of  providing  initial  values  for 
on-line  data  processing.  A  good  strategy  is  characterized  by  the  saving  of  process  time  and 
by  a  minimum  of  user  interaction,  such  as  for  for  choosing  the  method  to  be  used  or,  judg¬ 
ing  the  results  of  the  computations.  So,  for  the  rare  case  of  having  known  exterior  orienta¬ 
tion  elements  of  all  camera  stations  only  spatial  intersections  for  determining  object  point 
coordinates  are  used.  But  if  only  some  control  points  and  no  known  exterior  orientation  el¬ 
ements  are  available  (normal  case)  in  the  photogrammetric  block,  the  triangulation  has  to 
start  in  a  local  coordinate  system  by  orienting  the  first  two  images  relatively.  The  exterior 
orientation  elements  of  the  following  consecutive  images  in  the  strip  can  be  determined  by 
space  resection  if  sufficient  tie  points  can  be  nreasured.  The  transformation  into  the  global 
net  can  be  achieved  when  sufficient  control  point  information  is  available  (i.e.  for  the  three- 
dimensional  similarity  transformation).  In  an  idea)  case,  as  mentioned,  when  sufficient  con¬ 
trol  points  and  known  exterior  orientation  elements  are  given,  only  spatial  intersections 
have  to  use  for  confuting  object  point  coordinates  of  additional  nati^  points.  If  no  infor¬ 
mation  of  a  global  net  is  available  (bad  case),  the  triangulation  has  to  be  processed  in  a  lo¬ 
cal  datum.  Herein,  the  same  strategy,  as  used  for  a  global  datum,  for  providing  initial  values 
will  be  used  except  that  a  transformation  is  not  performed. 

OLTRIS  computes  initial  values  automatically  after  each  measurement,  and  before  updat¬ 
ing  the  normal  equation  system  with  any  new  observations,  if  sufficient  irtformation  for 
computing  initial  values  is  present  On  the  other  hand,  the  necessary  initial  values  of  all 
new  observations  will  be  computed,  if  the  operator  selects  the  option  Display  at  Update  on 
the  on-line  triangulation  processing  menu  (see  Figure  A  in  Appendix  A).  If  the  computa¬ 
tion  of  initial  values  is  not  possible  due  to  insufficient  information,  the  observations  to  be 
included  receive  the  status  waiting  at  not  included. 

The  fast  and  precise  computation  of  initial  values  improves  the  performance  of  an  on-line 
triangulation  and  data  processing  system  like  OLTRIS. 

5.2.  Adjustment 

Since  near  real-time  responses  are  crucial  factors  in  on-line  triangulation,  fast  sequential 
adjustment  algorithm  is  essential.  As  each  estimation  algorithm  can  be  reformulated  into  a 
sequential  mode,  a  great  variety  of  sequential  techniques  to  be  considered  arises.  The  im¬ 
portant  sequential  algorithms  and  also  a  brief  comparison  thereofis  presented  in  section 
5.2.1.  The  conclusion  was  to  implement  the  Givens  Transformations  (GT)  algorithm  for 
sequential  adjustment.  Because  of  avoiding  the  time  consuming  inversion  of  a  normal 
equation  system,  the  Cholesky  factorization  and  back  substitution  for  die  relinearizing  and 
update  of  the  parameter  vector  was  inqilemented,  for  the  simultaneous  least  squares  adjust¬ 
ment  case. 

5  J.l.  Simultaneous  adjustment  by  Choledty  factorization  and  back  substitution 

In  order  to  adjust  the  observations  simulttuieously  the  Cholesky  algorithm  was  implement¬ 
ed  in  OLTRIS.  This  algorithm  relinearizes  the  normal  equation  syston  and  updates  the  pa¬ 
rameter  vector. 

In  the  bundle  adjustment  it  is  assumed  that  the  linearized  observation  equation  system  is 
V  =  Ax-l;P  (4) 


KeRUn,  Holm,  Oruen 
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with: 


V . vector  of  residuals, 

A . design  matrix, 

X . vector  of  unknown  parameters, 

1 . vector  of  observadons  and 

P . weight  matrix. 


This  equation  system  is  solved  by  the  least  squares  method,  by  building  and  solving  the 
normal  equation  system  (5), 

a'^PAx  =  a'^PI  (5) 


or 

Nx  =  n  (6) 

with: 

N . normal  equation  matrix  and 

n . vector  of  the  right-hand-side. 


The  explicit  solution  to  the  unknown  parameters  may  be  written  as 


X  =  N"‘n  =  (A'^PA)"‘a‘^P1.  (7) 

Instead  of  direct  inversion,  the  Qiolesky’s  algorithm  decomposes  the  normal  equation  ma¬ 
trix  N  into  a  lower  triangular  matrix  C  and  an  upper  triangular  matrix 

N  =  CC'*'.  (8) 

Further,  the  vector  of  the  right-hand-side,  n  can  be  decomposed  as 

n  =  Cc  (9) 

The  reduction  (factorization)  of  the  normal  equation  matrix  N  yields  the  upper  triangular 
matrix  and  vector  c.  The  unknown  vector  x  is  solved  by  back  substitution  using 

C'^x  =  c.  (10) 
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Figure  S  displays  the  tnain  stq>s  of  Chole^  facttnization  and  back  substitution  graphical¬ 
ly. 


update  of  parameter  vector  by  Cholesky  algorithm 


5  J.2.  Sequential  atUustment  by  Givois  Transformations 

A  number  of  algorithms  for  sequential  adjustment  have  been  suggested.  Mikhail  and  Helm- 
ering  (1973)  present  methods  are  based  on  updating  the  inverted  normals  directly.  This 
method  is  also  known  as  “stationary  Kalman  filtering”  and  has  been  adopted  for  OLT  by 
several  authors  [Helmering,  1977;  Kratky,  1982;  Dowideit,  1980, 1982].  Chuen  (1982)  rec¬ 
ommended  the  lUangular  Factor  Updt^  (TFU),  which  updates  the  fimtorized  normals 
based  on  Gauss  or  Cholesky  decomposition.  The  procedure  of  TFU  is  also  described  by 
Wyatt  (1982).  He  compared  the  Kalrnan-fcam  to  the  TFU  algorithm.  The  ctmclusion  of  the 
comparison  was  that  Ae  ITU  algorithm  is  superior  to  the  Kalman-form  algoritiim  both 
with  respect  to  time  consumption  and  storage  requiremmts.  The  use  of  the  Givens  TYans- 
formations  in  photogrartunettic  applications  has  been  suggested  by  Blais  (1983). 

The  relative  performance  of  these  sequential  algorithms  can  best  be  judged  by  practical 
tests  using  software  which  is  as  close  as  possiUe  to  a  real  application.  Such  tests  were  car¬ 
ried  out  in  on-line  phototriangulation  by  Runge  (1987)  and  Holm  (1989c).  They  describe 
the  use  GT  in  on-line  photottiangulation  and  compare  this  approach  to  ibe  IHl  algo¬ 
rithm.  Range’s  results  of  die  conqiarison  with  the  TI^  method  showed  that  only  4)%  of 
the  time  is  needed  for  the  update  ctf  dw  reduced  normal  equations  when  using  GT.  In 
Holm’s  test,  the  GT  algorithm  qipears  to  be  up  to  four  times  faster  dian  the  TFU  for  iqxlat- 
ing  die  fiunotized  normal  equation  system. 

The  perfiotmance  of  GT  in  sequential  estimation  will  be  given  in  die  ftdlowing.  For  thme 
investigations  the  testfield  blo^  Echallens,  widi  a  total  of  49  jdMtt^iqdi  and  231  otiject 
points,  was  used.  The  conqnitation  was  made  on  a  SUN  3  computer  udng  a  68020  CPU 
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with  floating  point  accelerattx'.  (Note  that  the  computation  on  a  SUN  SPARCstation  l-i-  is 
30%  faster  than  on  the  SUN  3).  Figure  7  illustrates  the  results  of  the  conqiutations  of  the 
soluti(»  vectOT  in  a  sequential  mode  when  building  up  the  block  from  2  to  49  photognqihs. 
Herein,  groups  of  measurements  have  been  introduce  at  different  stages  ranging  from  34 
to  SO  image  points  per  update.  The  abscissa  gives  the  stage  of  the  sequential  process  (num¬ 
ber  of  photographs)  while  the  ordinate  gives  the  computing  time  per  new  image  point  in 
seconds. 


.•••.  / 


Hgure  6:  Block  Echaliens,  Computing  times  (solution  vector)  per  additional  image 
point  measurement  when  introducing  groups  of  measurements,  using 
Givens  transformations  algorithm  on  a  SUN  3  computer  upgrated  widt  a 
68020  CPU  with  floating  point  acceleratOT 

Further  information  about  this  test  of  the  performance  of  GT  can  be  found  in  the  second  In¬ 
terim  Report  [Holm,  1989b]. 

An  introduction  of  GT  describes  Blais  (1983): 

“The  Givens  Transformations  provide  a  direct  method  for  solving  linear  least-squares 
problems  without  forming  the  turrmal  equations.  One  most  important  feature  cf  this 
approach  for  recursive  applications  is  that  the  solution  can  be  obtained  at  any  stage 
(tfthe  processing  of  the  observational  irformation  by  simply  carrying  out  a  back  sub¬ 
stitution  for  the  unknown  parameter  into  an  upper  triangular  system  of  equations." 

Practically,  the  Givens  method  is  an  orthogonal  transformation  technique  based  on  the  use 
of  two-dimenskmal  rotations  to  eliminate  matrix  elements.  This  qiproach  is  a  fecial  case 
of  the  ^  decomposition  and  is  comparable  widi  the  CStoledcy  medred  as  described  in  sec¬ 
tion  5.2.1.  The  Givens  Iksnsflirmations  process  one  row  of  design  matrix  A  at  a  tinK 
and  are  used  for  the  sequential  addition  or  deletion  of  observations  in  an  interactive  envi- 
ronmenL  As  mentioned,  the  solution  vector  is  available  at  any  stage  by  sinqrly  performing 
a  back  substitution  in  an  upper  triangular  equatkra  sysieoa.  The  aflxemrationed  onbogonal 
decomposition  of  tiie  design  matrix  A  yields  Q  A  ■>  R.  Herrin,  the  upper  iriangnlar  mutix 
R  is  identical  widi  die  result  ot  the  Chriesky  decooqwrition  of  the  normal  equation  qrstem 
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ATPA.  When  A  is  a  m  x  n  matrix,  Q  can  be  interpreted  as  a  sequence  of  m  x  m  orthogonal 
matrices,  each  containing  transformation  coefiScients  of  one  single  GT. 

Range  (1987)  illustrates  the  functioning  of  the  GT  algmithm  using  an  exanqtle  of  a  3  x  3 
photo  block.  In  diis  example,  how  an  additional  measuremrat  of  a  new  image  point  yields  a 
pair  of  new  tow  vectors  containing  the  partial  derivations  of  the  new  measuranent  aixl  the 
respective  constant  values  is  denmnstnued.  This  new  observation  -  two  image  coordinates 
in  two  images  -  builds  the  four  rows  of  the  new  row  vector.  Individual  transformations  of 
these  four  tows  into  zero  cause  alterations  in  components  of  the  R  matrix.  Furthermore,  a 
transformation  of  the  respective  constant  values  yields  the  root  of  the  residual  sum  of 
squares.  After  the  transformations  the  whole  new  row  vector  contains  zeros  and  can  be 
fiUed  by  the  next  observation.  Sintilarly,  false  observations  can  be  removed  from  the  re¬ 
duced  normals,  the  constant  vectn,  and  the  square  toot  of  the  residual  sum  of  squares  by 
GT.  When  the  transfcnmation  coefScients  have  been  computed  with  the  diagonal  element 
of  j-th  row  of  matrix  R  and  of  the  design  matrix,  the  elements  of  matrix  R  and  A  will  be 
transformed  according  to  the  removed  observation.  GeneraUy,  after  updating  or  “downdat¬ 
ing”  the  reduced  normals  with  GT,  the  solution  vectrn:  can  be  comput^  with  hack  substitu¬ 
tion. 

Detailed  descriptions  and  also  mathematical  formulations  of  the  sequential  performaixre  of 
GT  are  given  in  Gruen  (198S),  Range  (1987),  Holm  (1989c)  and  Edmundson  (1991). 

5,3.  Quality  control 

The  major  goal  of  the  on-line  triangulation  techruque  is  to  control  tiie  quality  of  the  ob¬ 
served  data  at  an  early  stage  of  data  processing.  Tl^  permits,  in  the  presence  of  detected 
errors,  direct  re-measurement  of  obs^atkms  and  provides  a  cleaned  rfatji  set  for  the  final 
execution  of  the  simultaneous  block  adjustment  In  this  context  quality  control  consists  of 
two  tasks:  (a)  detection  of  blundns,  using  Baatda’s  data  snooping  techruque  [Baarda, 
1968].  (b)  compensation  of  the  systematic  errors,  using  additional  parameters  for  self  cali¬ 
bration. 

5  Blunder  detection  and  locatitm 

^thin  the  least  squares  estimation  qrproach,  Baarda’s  data  snooping  technique  has  proven 
to  be  an  effective  method  for  blunder  Election.  This  method  has  been  suggested  for  on-UiK 
triangulation  by  Foersmer  (1979),  Kratky  (1980b),  Moloiaar  (1981c),  Dowideit  (1982)  and 
Gruen  (1982).  The  m^Kxl  requires  only  the  diagonal  elements  of  the  matrix  (weight 
coefficient  matrix  of  the  residuals) 

Q,»  =  Qii-AQ„A'f  (11) 


with: 

(2d . weight  coefficient  matrix  of  the  observations, 

(2xit . weight  coefficient  matrix  cf  the  unknowns  and 

A . design  matrix. 


A  ranarkable  computational  speed-up  of  the  datt  snooping  techruque  can  be  achieved  with 
the  method  of  “unit  observation  vector^  [Gruen,  1982],  if  cmly  a  few  observations  have  to 
be  tested  at  any  given  time  and  only  the  related  diagonal  efemmts  of  the  <2w  nutiix  are  re- 
quiied.  This  is  likely  to  be  die  case  in  on-line  trianguluion,  where  observations  acquired  at 
earlier  stages  the  seqtrential  process  have  already  passed  qnalhy  oontred  tests. 
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Baaida  suggested  the  test  criterion  Wj 

|Vi| 


with: 


V. 

where: 

Wj^ . standardized  residual, 

o„ . estimated  variance  of  unit  weight  and 

Vj . residual, 

. estimated  variance  cd  residuals, 

q,  . i*  diagonal  element  of  Qw  -matrix. 


(12) 


Under  the  null-hypothesis  :  E(vi)  =  0,  W|  is  noimaUy  distrilnited,  so  that  critical  values 
from  the  normal  table  can  be  ^sed  for  testing.  In  practice,  however,  the  expected  value  of 
the  vtuia^  of  unit  weight  {t^is  not  available,  lltus,  an  estimate  of  die  variance  of  unit 
weight  dgfiom  the  simultaneous  adjustment  or  from  any  sequential  stage  must  be  used  for 
testing  standardized  residuals. 

Under  Hoi ,  wj  now  has  a  student-t  distribution  and  again  critical  values  for  testing  are  easi¬ 
ly  available  from  the  respective  tabtes.  For  exairqtle,  at  a  a  0.001  ctmfidence  level,  (giv¬ 
ing  a  probability  of  99.9%),  and  assuming  an  infinite  number  of  observations,  the  critical 
test  value  is  c  =  3  JI91  fix’  the  detection  of  a  gross  etror. 

Note  that  it  is  assumed  that  only  <me  blunder  is  present  in  tite  photogrammetric  system.  But 
practical  experience  shows,  however,  that  such  an  assimqttion  is  not  realistic.  In  addition  to 
this  proUein,  correlations  between  the  residuals  cause  the  effects  of  blunders  to  be  either 
conqiensated  or  accumulated,  malting  blunder  detection  and  location  a  much  more  compli¬ 
cated  affair. 

To  perfixm  reliable  Uunder  detection,  it  is  necessary  to  add  more  than  the  traditional  two 
photos  a  stereo  model  into  the  normal  equation  system.  For  higher  redundancy  and  better 
reliability,  a  larger  number  of  nys  to  each  object  prnnt  is  required.  Ihus,  more  dian  two 
photos  should  be  treated  as  an  unit,  i.e.  congaited  as  a  snbblo^  For  an  object  pctint,  mea¬ 
sured  with  two  rays,  a  Uunder  in  the  x-image  coordinate  is  not  deiectaUe  and  in  y-image 
coorditute  not  locatable.  A  Uunder  in  the  x-measuremeat  of  a  duee  ray  object  pram  is  also 
not  locatable.  Consequendy,  reliable  and  efficient  Uunder  detection  location  requires 
more  than  three  rays  to  an  objea  pUnt  According  to  Omen  (1982),  if  the  overlap  is  60%,  a 
subUock  of  3  X  3  photos  gives  optimal  local  reliability  in  aerotiiangulation. 

PiacticaUy,  the  critical  aspect  from  a  hardware  perspective  is  die  speed^perfixmance  of  the 
conqwtadon  of  a  test  criterion  W}  fix  each  observation.  Theiefixe  investigations  into  the 
computation  of  the  test  criteria  haire  been  dcxie  under  die  same  conditioos  as  mentioned  in 
5.2. 1 .  Rgme  8  diows  avetife  ootqmting  times  fix  Wi  per  image  point  with  respect  to  Uodc 
sine.  For  this  purpose,  the  test  Uock  Ecfaallens  was  as  a  cotqdere  Uock  (7  X  7  photos) 

md  two  snbbkick  nudes  (3  X  3  and  S  X  5  photos).  In  eadi  case,  a  wieqr  of  image  point  ob¬ 
servations,  tanghtg  Cram  2  w  1736  observations,  was  tested  nring  dw  data  snooping  test 
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Figure  7:  Block  Echallens.  Computing  times  for  W}  per  image  point  with  respea 
to  block  size  on  a  SUN  3  cmi^uter  (68020  CPU  wi±  floating  point  ac¬ 
celerator) 


criteria  W|.  It  can  be  seen  that  the  increase  of  computing  time  is  nonlinear  with  respect  to  an 
increase  in  the  block  size. 


5J^  STKcmatk  error  compensation 

The  presence  of  systematic  image  eirars  in  photogrammetric  blocks  is  a  serious  problem. 
Systematic  errors  can  be  caused  by  the  folkrwing  effects: 

•  Atmospheric  refraction 

•  Geometric  distortion  by  imperfections  of  lenses 

•  Change  of  camera  constant  and  princ^le  point  by  variation  of  die  temperature  and 
pressure  during  photognqifaing 

•  Inqnecise  fiducial  marks 

•  Vacuum  plate  unflatness 

•  Him  and  stratum  deformation 

•  Blur  effected  by  image  motion 

Different  qiproaches  can  be  enqiloyed  to  eliminate  andAv  conqiensaie  these  systematic  er¬ 
rors.  The  most  efficient  method  is  the  simultaneous  conqiensation  of  the  systematic  errors 
expanding  the  estimation  model  with  additional  parameteo.  This  procedure  is  known  u 
“self  calibruion’*.  Self  caUbrmioo  has  proved  its  potential  in  off-line  triangnlation.  If  on¬ 
line  tringnlation  is  petfocmed  a  high  accuracy  level,  self  calflntion  is  also  indiapens- 
dde. 
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TVvo  sets  cd  parameters  are  suggested  for  the  application  of  aerMriangulation; 

•  12-parameta'  set  of  Ebner  (1976).  an  “orthogonal  moder  referring  to  3  x  3  image 
pcwt  distribution 

•  44-parameter  set  of  Giuen  (1978),  an  “mthogonal  model”  referring  to  5  x  5  image 
point  distributicm 

The  parartxter  set  of  Ebner  has  been  used  for  cartographic  and  upognphic  aerotriangula- 
tion,  while  Gnien’s  parameter  set  is  relevant  for  apphcaticms  in  cadastral  surveying,  etc. 

In  the  initial  stages  of  a  sequential  OLT  system  the  probability  of  carrying  non-determin- 
able  additioiud  parameters  is  certainly  higher  than  in  an  off-line  ^tem.  Non-detenninable 
additional  parameters  can  significantly  weaken  the  system  and,  in  extreme  cases,  lead  to 
near-singuiirity.  In  mder  to  avmd  running  into  too  many  non-determiiuible  additional  pa¬ 
rameters  and  spending  too  much  rime  with  the  clearance  procedure,  it  is  advisable  to  oper¬ 
ate  an  OLT  with  a  fairly  small  set  of  additional  parameters,  considering  only  the  most 
important  ones  which  ate  known  to  be  determinable  in  those  relevant  arrangements.  A 
computational  fast  and  reliable  procedure  for  the  handling  of  additional  parameters  (check¬ 
ing  on  detenninability  and  possibly  significance)  has  still  to  be  adopi^  for  OLT.  Such  a 
procedure  has  been  suggested  by  Gruoi  (1978a,  1983a),  Foerstner  (1981)  and  Jacobsen 
(1982).  The  summarized  mathen^cal  fixmulation  of  different  sets  of  additional  parame¬ 
ters  has  been  published  in  Kilpell  (1980). 

Although  preparatory  studies  for  self  calilsatimi  have  been  made  on  this  project,  due  to  the 
redirectioa  of  the  other  project  wmk  and  the  resulting  time  constraints,  its  implementation 
into  OLTRIS  could  not  be  ctxnpleted. 

6.  Software  modules  of  OLTRIS 

The  OLT  system  is  decomposed  into  several  modules.  Hgure  8  gives  an  overview  of  the 
main  modules  with  their  interrelations  and  dependencies.  The  shape  of  die  boxes  in  the  di¬ 
agrams  denotes  whether  the  modules  contain  mainly  functions  (rectangles)  ex’  are  primarily 
used  for  data  storage  (rounded).  The  following  modules  contain  the  administration  of  Ol^ 
TRIS  (top  level): 

•  USER_INT/XVIEW;  User  interface  handles  the  communicaticHi  between  the  user  and 
the  systertL  The  graphical  user  interface  is  based  on  OpenWindows/OPEN  LCXJK. 
Rfithin  this  user  inioftce  commands  ate  to  be  given  by  menus,  buttons,  and  text  in 
windows. 

•  FLOW_CXH4:  This  administration  module  contains  the  flow  control  of  die  systnn  for 
the  non  grqihical  user  interface.  The  flow  control  for  die  program  administrated  by  the 
gnqihical  user  interface  is  handled  in  die  module  XVIEW. 

The  cote  functions  are  contained  in  die  second  level  of  modules. 

•  INPUT  The  functions  of  diis  module  will  control  the  data  input  from  outside  of  the 
system,  osganize  and  store  it  in  internal  data  structures.  Data  may  come  from  the  user 
dialogue,  measurements  (in  the  final  verskm)  or  from  files. 

•  INITIAL:  The  computation  of  initial  values  for  all  unknown  parameters  is  necessary 
becaoK  of  the  non-linear  equation  systera  Diffineot  frmcDons  have  been  inqihnnented 
for  computing  the  initial  values  of  object  points  and  exterior  oifentation  parameters. 

•  ADJUST  This  adjustment  module  contains  frmetions  for  die  sjnmltaneous  adjustment 
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Hgure  8:  Software  modules  of  (be  OLTRIS. 


(Cholesky  factorization  and  back  substitution)  and  sequential  estimation  (Givens 
transformations). 

•  QUAL_CX)N:  The  quality  control  module  contains  fimctions  for  computing  some  gen¬ 
eral  statistics,  summarizing  the  results  as  a  whole,  performing  the  blunder  detection 
with  data-$noq)ing.  It  shoi^  be  expanded  with  fimctions  for  checking  determinability 
and  significance  of  additional  parameters. 

•  REPORT:  Functions  for  the  output  of  results  and  data  as  hard  copies  or  on  secondary 
storage  for  further  use  are  found  in  this  module. 

The  following  two  modules  contain  a  library  of  conunon  functions: 

•  FUNC_LIB:  This  module  contains  a  coUecdmi  (library)  for  ctmimon  functions  used 
by  dififerent  modules.  In  particular,  the  functions  are  for  computing  mathematical  ap¬ 
plications,  e.g.  rotation  matrix,  design  matrix  coefficients,  and  for  administrating  com¬ 
mon  a^ilications,  e.g.  input  of  values,  file  handling  etc. 

•  MEH_ALLOC:  Functions  for  memory  allocatioo  and  reallocation  of  specific  items  of 
the  data  structures  defined  in  the  system  are  collected  in  fliis  module. 

The  following  modules  contain  mainly  data: 

•  CONTrLPAR:  This  nxxhile  contains  cmitiol  flow  parameters,  as  well  as  critical  val¬ 
ues  and  other  general  parameters  needed  some  conqnitations. 

>  Mffl^LJ^ATA:  Ctata  for  entities  in  the  madiematical  modd,  e.g.  cameras,  images, 
object  points,  alio  a  priori  and  conyuted  parameters,  obsCTvations  and  dedgn  matrix. 


Pirn 
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are  included  in  this  nxxlule. 

•  EQ_SYST:  This  module  includes  data  and  functions  fm  handling  the  nmnial  equaticm 
system  and  related  data. 

Detailed  description  about  the  modules  and  its  containing  functions  is  given  in  the  software 
documentation  (Holm,  Kersten,  1991]. 

7.  The  program  OLTRIS 

OLTRIS  has  been  designed  as  an  on-line  data  processing  system  for  measuring  image  point 
coordinates  and  adjusting  the  observations  sequentially  and/or  simultaneously.  Additional¬ 
ly,  object  point  cocndinates  and  also  exterior  orientation  elements  can  be  included  as  obso*- 
vations.  Image  point  coordinates  can  be  measured  in  mono  or  stereo  mode.  Furthermore, 
the  system  provides  blunder  detection  with  data  snooping  after  the  adjustment.  Finally,  out¬ 
put  fiinctionallity  (e.g.  displaying  of  information,  plotting  and  saving  of  data  and  of  results) 
is  provided  in  OLTRIS. 

7.1.  The  user  interface 

The  program  is  to  be  processed  by  menus,  buttons,  and  text  in  windows  of  a  graphical  user 
interface.  The  user  interface  is  based  on  OpenWindows/OPEN  LOOK.  A  prevoius  key¬ 
board/terminal  version  of  the  program  results  from  the  first  implementations  of  the  source 
code,  and  was  used  mainly  for  testing  purposes.  An  example  for  keyboard  dialog  of  the 
program  was  included  as  an  tqipendix  in  the  6th  Interim  Report  The  graphical  user  inter¬ 
face  has  been  developed  with  the  OpenWindows  Developer’s  Guide  (Devguide)  2.0,  with 
more  additions  under  version  3.0.  The  following  discussion  concentrates  on  the  use  of  the 
program  using  this  graphical  interface. 

The  data  processing  in  OLTRIS  can  be  separated  in  four  phases:  initialization,  measure¬ 
ment  and  data  processing,  quality  control,  and  documentation.  The  last  three  phases  cerate 
in  parallel. 

7.2.  Initialization 

The  first  phase  (initialization)  ctmtains  the  preparation  of  OLTRIS  for  on-line  data  process¬ 
ing,  which  means,  in  general,  the  input  of  data  and  the  setting  of  control  parameters.  In  par¬ 
ticular,  OLTRIS  provides  fm  the  input  of  control  parameters  (e.g.  iteration  stop  criterion, 
maximum  numbCT  of  iterations  etc.),  camera  parameters  (e.g  camera  calibration  data  and 
related  informations),  the  measured  inuige  coordinate  oteervadons,  object  point  comdi- 
nates  and  exterior  orientation  elements,  using  the  last  two  mentioned  input  options  as  ob¬ 
servations  or  initial  values.  The  options  for  initialization  -  its  panel  can  be  started  frtnn  the 
on-line  triangulation  processing  panel  (see  Figure  A  ctf  Append  A)  -  are  presented  in  Fig¬ 
ure  B  of  Appendix  B.  In  initiali^on,  the  observations  of  image  pmnt  coordinates,  object 
point  coontoates  and  exterkH*  orientation  elements  are  read  from  ^  and  ate  placed  in  tem¬ 
porary  data  structures  of  die  system  (intenul  buffer).  During  die  measurement  procedure 
die  observation  data  are  taken  from  this  internal  bufi^  and  placed  in  the  data  structures  of 
the  normal  equation  system,  lb  date,  no  direct  input  from  a  digital  workstation  or  an  ana¬ 
lytical  plotter  is  not  yet  possible.  The  input  of  camera  calibration  data  (see  Hgure  D  in  Ap¬ 
pendix  Q  and  image  or  modd  data  (see  Rgure  E  in  ^ipendix  C)  is  possible  within  the 
initialization  menu  and  also  later  witl^  the  input  option  of  the  on-line  ttiangnlation  ]m>- 
cessing  pmeL  Rgure  9  gives  an  overview  of  die  data  flow  in  the  system.  The  Idk  side  of  the 
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Figure  9:  OLTRIS  and  its  data  flow 


figure  shows  the  input  possibilities  into  the  system,  while  on  the  right  the  output  of  results 
and  data  is  presented.  Itetween  processing  (as  far  as  the  data  are  processed),  die  results  and 
also  the  processed  data  can  be  displayed  and  plotted  cm  screen,  reported  on  paper  and/or 
stored  on  disk.  Automatic  ccmtiol  be  achieved  by  setting  control  flow  parameters  -  Fig¬ 
ure  C  in  Appendix  B  illustrates  its  panel  -  and  the  user  interactions  lead  the  processing  of 
the  system. 


7J.  Measurement  and  data  processing 


The  second  phase  of  data  processing  in  OLTRIS  is  for  image  measurement  and  data  pro- 
ces^g.  When  starting  up  the  program,  the  on-line  triangulation  processing  panel  (main 
panel)  iqrpears  (see  Figure  A  in  Appendix  A)  on  screen,  from  which  all  processes  can  be 
managed.  After  initializatiem  of  parameters  arxl  data  input  as  mentioned  befene.  the  mea¬ 
surement  jnocedures  can  be  started  frtnn  the  main  menu.  The  following  possibilities  for 
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Hgnre  10;  Measuremoit  options  of  (XITRIS 

measuring  image  point  coofdinaies  as  depicted  in  Rgure  10  can  be  sdected  by  the  operator. 
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When  begmiiing  a  new  triangulation  it  is  advisable  firstly  to  measure  in  the  start  image  (or 
model)  all  available  control  points  in  order  to  fix  the  datum  and  secondly  to  continue  with 
“next  observation”  or  "next  and  more”.  Continuing  triangulation  by  using  the  option  “tie 
points”  all  following  images  (or  models)  will  be  connected  in  the  block  sequentially.  The 
measurement  options  “One  specific”  and  “Manually”  provide  an  interaction  of  the  cqierator 
for  measuring  cotain  points.  The  measured  points  will  be  plotted  on-line  in  a  window  (see 
Figure  I  in  Appendix  E)  to  show  the  coniributicm  of  points  in  the  images.  After  each  mea¬ 
surement  the  normal  equation  system  is  updated  sequentially  with  the  new  observations,  if 
initial  values  for  the  respective  object  point  and  exterior  orientation  parameters  of  the  new¬ 
ly  measured  image  are  available.  Otherwise  the  iiutial  values  will  be  ccniqiuted  when  dis¬ 
playing  data  or  telinearizing  and  updating  the  equation  system  by  pressing  the  button 
Display  or  Update.  The  functional  co-operaticm  of  the  sequential  and  simultaneous  relin- 


Figure  11:  Relinearizing  and  updating  of  the  normals 


by  Givens  ‘Ihuisformations  and  (!holesky 

earization  and  update  of  the  normal  equation  syston,  udien  introducing  new  (additional) 
observations,  is  illustrated  in  Hgure  11.  The  strategy  and  thefunctitms  tocranputing  initial 
values  are  desoibed  in  section  5.1.  The  adjustment  procedures  are  introduce  in  section 
5.2.  To  control  the  adjustment  the  operator  has  to  choose  die  iteration  criterion  (avoage 
changes  of  the  rotation  angles  or  of  the  coordinates,  or  the  average  changes  of  v'^pv),  the  it- 
oation  stop  criterion  and  die  maximum  number  iteiatkms  (in  initialization). 
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7.4.  Quality  control 

During  measurement  and  data  processing  OLTRIS  provides  a  quality  control  of  the  obser¬ 
vations  with  data  snooping.  The  goal  of  the  quality  control  is  the  detection  and  location  of 
blunders  in  the  observed  data.  To  this  end.  in  addition  to  image  pmnt  coordinate  observa¬ 
tions,  all  available  observations  of  exterior  orientation  elements  and  object  points  will  be 
tested  by  the  data  snooping  method.  An  example  of  a  detected  blunder  is  shown  in  Figure  F 
(Appendix  D).  For  cleaning  the  data  set,  a  detected  blunder  can  be  excluded  from  the  ad¬ 
justment  It  is  recommended  to  remove  the  detected  observation  finally  from  the  internal 
buffer  only  after  performing  a  simultaneous  adjustment  in  order  to  the  see  the  effect  of  its 
exclusion. 

IS.  Documentation 

During  and  after  measurement  the  operator  may  want  to  see  up-to-date  results  and  data 
records.  Therefore,  the  third  phase  (documentation)  cannot  be  separated  from  the  measure¬ 
ment  and  data  process  phase.  The  documentation  of  results  and  processed  data  is  separated 
in  Display,  its  options  are  depicted  in  Ingure  12,  and  in  Report,  its  options  am  Save,  Print 
and  Plot.  The  panels  of  these  report  options  are  shown  in  Appendix  D,  Figures  F,G,H. 


Figure  12;  Display  optitms  of  OLTRIS 


7.6.  Results 

To  illustrate  the  performance  and  functionary  of  OLTRIS,  a  step-by-step  triangulation  of  a 
photogrammetric  block  is  demonstrated  stepwise.  The  data  of  the  photogrammetric  block 
Simplon,  a  mountainous  area,  has  been  chosen  far  the  demonstration.  In  order  to  reduce  dte 
quantity  of  shown  data  and  results  m  a  minimum  in  the  report,  only  die  start  and  the  end  of 
the  triangulation  will  be  presented.  The  chosen  block  was  flown  in  four  strips  with  60% 
overlap  in  strip  and  cross-strip  direction.  The  processing  of  the  block  (4x6  images,  102 
object  points,  943  observations,  450  unknowns  and  493  redundant  observaticnis)  was  exe¬ 
cuted  on  a  Sun  SpaicStatitm  1.  The  first  seven  exterior  orientatiem  elements  are  introduced 
as  observations  thereby  fixing  the  datum  of  the  photogrammetric  block.  After  introducing 
the  two  first  images  and  measuring  6  points  in  each  image  (see  Hgure  J  in  Appendix  E),  at 
this  stage,  the  adjustment  can  be  computed  simultaneously  for  efficiency  and  stability  rea¬ 
sons.  Having  one  redundant  observation  (see  Rgure  K  in  Appoidix  E),  ^  sigma  a  posteri¬ 
ori  has  been  estimated  in  die  adjustment  as  11.7  microns,  compared  to  sigma  a  priiai 
d^ned  at  5.0  microns.  Adding  more  images  and  observations  increases  the  redundancy 
and  die  photogrammetric  Mock  becomes  more  stable.  The  perfcnmance  of  the  sequential 
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update  of  the  noraial  equation  system  when  introducing  new  observations  is  shown  in  Fig¬ 
ure  7,  but  processing  a  different  photogrammetric  test  block  on  a  Sun  3  workstation.  The 
result  of  the  final  adjustment  of  the  Simplon  block  is  presented  in  Figure  L,  in  Appendix  F. 
A  confirmation  of  the  quality  of  the  mathematical  model  is  the  value  of  sigma  a  posteriori, 
which,  as  seen,  does  not  differ  greatly  from  the  introduced  value  of  sigma  a  priori.  As  an 
example  of  the  results  the  solution  of  the  exterior  orientation  elements  is  depicted  in  Figure 
M  of  Appendix  F. 

In  summary,  a  first  version  of  the  system  has  been  implemented  on  the  Sun  workstation.  In 
the  future,  it  has  to  be  decided,  to  which  measurement  system  -  analytical  plotter  or  digital 
workstation  -  OLTRIS  should  be  connected.  Finally,  there  is  always  room  for  improvement 
of  a  very  complex  program  like  OLTRIS. 

8.  Conclusions 

To  the  authors’  knowledge,  no  fully  operational  OLT  system  has  been  realized  yet.  A  real¬ 
ization  of  a  practical  OLT  system  could  look  as  depict^  in  Figure  13.  The  left  side  of  the 
Figure  illustrates  the  input  into  the  system  of  direct  image  measurements  by  an  Analytical 
plotter  or  of  digital  images  by  an  image  acquisition  station  via  a  scanner  or  via  CCD-cam- 
eras.  The  data,  as  shown  on  the  right  side  of  Figure  13,  can  be  saved  on  disk  or  printed  as 
hard  copies.  Automatic  control  and  user  interactions  lead  the  processing  of  an  OLT  System. 

Theoretical  OLT  problems,  except  statistical  testing  in  sequential  mode,  have  been  solved, 
but  have  not  been  practically  combined  into  a  single  system.  The  continual  improvement  in 
computer  performance  has  caused  reduced  efforts  in  the  development  of  fast  processing  al¬ 
gorithms  for  OLT.  Furthermore,  investigations  into  instrument  development  are  concentrat¬ 
ing  on  the  improvement  of  automatic  nwasurcment  procedures  on  analytical  plotters - 
instead  of  implementing  fast  data  processing  algorithms,  including  quality  control  to  detect 
gross  errors  and  to  check  precision  and  reliability  of  the  network. 

The  authors  see  potential  in  further  development  and  use  of  the  OLT  algorithms  and  meth¬ 
ods  for  applications  in  robotics,  industrial  quality  control  and  autonomous  vehicle  naviga¬ 
tion.  Of  particular  interest  in  robotics  and  autonomous  vehicle  navigation  is,  for  example, 
the  task  of  determining  the  position  and  orientation  of  a  mobile  robot  on-line.  Here  on-line 
data  processing  could  also  play  an  important  role. 
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Figure  13:  Hardware  configuration  of  an  OLT  system 


10.  References 
10.1.  Internal  reports 

Kersten,  Th.,  Holm,  K.  R.,  1991.  Software  Documentation  (SD.OO  -  SD.03),  Vision  2.0, 
June  1991,  IGP/ETH-Zurich. 

Holm,  K.  R,  1989a.  1st  Interim  Repon  (Jan.-Feb.  89),  IGP/ETH-Zurich. 

Holm,  K.  R.,  1989b,  2ik1  Intoim  Report  (Mar.-Jun.  89),  IGP/ETH-Zurich. 

Kersten,  Th.,  1990a.  3rd  Interim  Report  (Jul.-Dec.  89),  IGP/ETH-Zurich. 

Kersten,  Th.,  1990b.  4th  Interim  Report  (Jan.-Jun.  90),  IGP/ETH-Zurich. 

Kersten,  Th.,  1991a.  3th  Interim  Report  (Jul.-Dec.  90),  IGP/ETH-Zurich. 

Kersten,  Th.,  1991b.  6th  Interim  Report  (Jan.-Jun.  91).  IGP/ETH-Zurich. 


Kenia,  Hdm,  Gnien 


Page  23 


On-line  point  posidMiing  with  single  fiame  camera  daB 


10  General  references 

Baaida,  W.,  1968.  A  testing  procedure  for  use  in  geodetic  networks.  Netherlands  Geod. 
Com.  Publications  on  Geodesy.  Vol.2,  No.S,  Delft 

Baarda,  W.,  1977.  Measures  of  the  accuracy  of  geodetic  networks.  lAG  Symposium  on  Op¬ 
timization  of  Design  an  Computation  of  Control  Networks,  Sopron,  pp.  419-436. 

Blais,  J.  A.  R.,1983.  Linear  least-squares  computations  using  Givens  transformation.  The 
Canadian  SurveyOT,  Vol.  37,  No.4,  pp.  225-233. 

Cogan,  L.,  1986.  Design  concepts  of  on-line  aerial  triangulation  software.  ASP/ACSM 
Convention,  Washington. 

Dorrer,  E.,  1978.  A  contribution  to  sequential  on-line  numerical  aerial  ttiangulation.  Pre¬ 
sented  Paper,  Syn^sium  of  the  Commission  in  of  ISP,  Moscow. 

Dowideit,  G.  R.,  1980.  On-line  bundle  adjustment  for  analytical  plotter.  Int  Archives  of 
Photogr.,  \bl.  XXm,  B3,  Commission  HI,  pp.  168-177. 

Dowideit,  G.  R.,1982.  Experience  with  improved  on-line-bundle-block-adjustment  on  the 
APC3.  Proceedings  of  the  Symposium  of  Com.  m  of  ISPRS  Helsinki,  June  7-11, 1982. 

Ebner,  H.,  1976.  Self  Calibrating  Block  Adjustment  Int  Arch.  Phot,  VoL  XXL  Part  3,  In¬ 
vited  Paper.  Com.  m,  ISP  Congress,  Helsinki. 

Edmundson,  K.  L.,  1991.  On-line  triangulation  of  sequential  stereo-pairs  using  Givens 
Transformations  without  square  roots.  Master  Thesis,  Department  of  Geodetic  Science 
and  Surveying,  The  Ohio  State  University. 

Foersmer,  W.,  1979.  Sequential  checking  of  model  coordinates  with  a  program  for  strip  ad¬ 
justment  Presented  Paper,  Symposium  on  On-Line  IViangulation,  ISP  WG  in/4  ,  Otta¬ 
wa. 

Gnien,  A.,  1978.  Progress  m  photogrammetric  point  determination  by  conqrenstion  of  sys¬ 
tematic  errors  and  detection  of  gross  errors.  Presented  paper  at  Symposium  of  the  Com¬ 
mission  m  of  ISP  Moscow,  1978. 

Gruen,  A.,  1982.  An  optimum  algorithm  for  on-line  triangulation.  Photogranunetric  Engi¬ 
neering  and  Remote  Sensing,  VoLSl,  No.4,  pp.419-436. 

Gruen,  A.,  1985.  Algorithmic  aspecu  in  on-line  oiangulation.  Paper  presented  to  the  Sym¬ 
posium  of  the  Commission  m  of  ISPRS  Helsinki,  June  7-11, 1982. 

Gruen,  A.,  1986a.  Photogrammetrische  Punktbestimmung  mit  der  Biindelmethode.  Institut 
fUr  GeodHsie  und  Photogrammetrie  an  der  ETH  Zurich,  Mitteilungen  Nr.  40. 

Helmering,  R.  J.,  1977.  A  general  sequential  algorithm  for  photogrammetric  on-line  pro¬ 
cessing.  Photogramm.  Eng.  and  RenxKe  Sensing,  M>L  43,  pp.  469-474. 

Holm,  K.  R.,  1989c.  Test  of  algorithm  for  sequential  adjustment  in  on-line  triangulation. 
Photogrammetria(PRS),  43,  pp.  143-156. 

Kratky,  V.,  1979.  On-line  analytical  triangulation.  The  Canadian  Surveyor,  Vol.  33,  No.2. 

Kratky,  V.,  1982.  Present  status  of  on-line  analytical  ttiangulation.  Photogtammetria(PRS), 
38,  pp.1-12. 

MikhaiL  E.  M.,  Helmering,  R.  J.,  1973.  Recursive  methods  in  photogrammetric  datareduc- 
ti<m.  Photogramm.  Eng.,  39(9),  pp.  983-989. 

Runge,  A.,  1987.  The  use  of  Givens  Transformation  in  on-line  photottiangulation.  Pro¬ 
ceedings  of  Intetcommission  Conference  m  Fast  Processing  of  Photogrammetric  Data, 
Interlaken,  Switzeiland,  June  2-4, 1987,  pp.  179-203. 

l^^tt,  A.  R,  1982.  On-line  photograminetric  triangulation  -  An  algwithmic  apfvoach. 
Master  Thesis,  Department  of  Geodetic  Science  and  Surveying,  The  Ohio  State  Univer¬ 
sity. 


On-line  point  positioaing  with  single  fiaune  camera  data 


Appendix  A 


A|>pendixB 


On-line  pcmt  positioniiig  with  single  bame  camera  data 


RgureC:  Subpanel  of  the  on-line  triangulation  initialization  panel 


Figure  D:  Input  of  camera  calibration  data  (Subpanel  of  the 
on-line  triangulation  initialization  panel) 


HguieE:  Input  of  image  or  model  data  (Subpanel  of  die 
on-line  triangnlatkm  initialintion  panel) 
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On-line  poiiitpositioiuiig  with  single  fame  c«ineni«<««» 


RgureF:  Display  for  quality  coitrol;  detected  blunder  in  point  2021  ofimage210 
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HguieQ:  Panel  for  printing  data  and  results 


HgureH:  Panel  for  saving  data  and  results 


HgureK:  Diqdiy  triangolitkm  stttu  after  the  fint  modet 
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SOFTWARE  DOCUMENTATION  SCHEME 


1.  Scope  of  the  documentation 

This  software  documcnatioo  is  based  on  the  recommendations  outlined  in  IEEE  1016.  A 
rigtxDus  "Software  Requirements  Specification"  as  described  in  IEEE  830  is  not  made  in  this 
pioject  The  requirements  for  the  scrftware  are  instead  embedded  in  this  descriptimi. 


2 .  Description  Organisation 

TF.EF.  1016  (Section  6.)  contains  four  different  methods  for  describing  the  design  of  software. 
Two  of  these  four  appn>aches,£>econv)iositu>n  and  Detailed  Design  have  been  chosen  for 
documenting  the  sof^are  of  OLTRIS. 

Decomposition:  Identification.  Type,  Purpose,  Function,  Submdinates  Dependecies, 
Resources. 

Detailed  Design:  Identification,  Function,  Interface,  Processing,  Data. 

The  characteristics  ai  the  design  description  are  defined  in  section  3. 

3.  Definitions  and  acronyms  used 

Design  Entity:  (See  IEEE  1016,  Sect  52) "...  an  element  (conqxment)  of  a  design  that  is 
stnicturally  and  functionally  (fistinct  fiom  odier  elements  and  that  is 
separately  named  and  re^nced".  Objective  of  decomposition  is  "to  divide 
the  system  into  separate  cooqxments  dhat  can  be  considered,  in^rlemented, 
changed,  and  tested  with  mininia!:  effect  on  other  entities". 

Entity  Attribute:  (See  IEEE  1016,  Sect  5.3) "...  named  characteristic  or  propeity  of  a  design 
enti^". 

The  individual  attributes  which  are  addressed  in  the  design  descriptions,  are 
listed  below,  with  shott(ened)  descriptions  (see  IEEE  1016, 5.3.1-5.3.10); 

1.  Identification  A  name;  unique;  duuacteiizing  the  nature  the  entity. 

2.  Type  Kind  of  entity.  Subprogram,  module,  procedure,  process,  data  store. 

In  this  documentation,  "module"  means  a  design  entity,  without  regarding 
how  it  will  be  inqrfemoiied,  while  "program  n^ule"  ideates  a  physical 
module  e.g.  in  a  C  program. 

3.  Purpose  A  descr^on  rfi  why  the  entity  exists  -  ...the  q>ecific  functional  and 

performance  requirenoents  for  which  diis  entity  was  created. 

4.  Function  A  statonent  oS  what  the  entity  does  - ...  the  transformation  applied  to  the 

inputs  to  ptt^uce  the  output  (Deodls  belong  in  9.  nocess). 

Fw  data  entity:  Type  of  information  stored  or  transferred. 

5.  Subordinates  Identification  of  all  entities  conqx>sing  this  one. 

6.  Dependencies  Descriptioos  of  dretelationsh^  with  other  entities.  Shall  identify  die 

"uses..."  or  "requires  the  presence  of ..."  relationship.  May  be  ducted  by 
structure  chans,  data  flow  diagrams,  transaction  diagrams. 

7.  Interfax  Descrqxionofhowoiherentitiesintetacis  with  this  entity.  Methods  and 

rules,  mcluding  e.g.  oonaiwnicating  through  parameirrs  and  direct  access  to 
intend  data  (Detute  in  IEEE  101^ 

S.Resources  Deacription  of  etemaits  that  are  exienwd  to  die  design.  Physical  devices, 

software  services  (e.g.  libtaries),  processing  lesouroes. 
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9.  Processing  Description  cf  the  rules  used  by  die  entity  to  adiieve  its  fiincdon.  Algoridun, 

contingencies.  A  refinement  of  the  Function  attribute.  Should  include: 

timing,  sequencing,  pierequisits  for  initiation,  priority  of  events . 

tennination  criteria,  handling  of  contingencies. 

10.  Data  Description  data  elements  internal  to  the  entity. 

Repie^tation,  fcmnat,  structure,  initial  values,  meaning  and  use, 
ac^table  values.  May  be  in  die  fonn  of  a  data  dictionaiy. 


4 .  References 

IFRF.  1016  IEEE  Std  1016-1987:  "IF.ER  Recommended  Practice  fOT  Software  Design 
Descriptions". 

IEEE  Std  830-1984:  "IEEE  Guide  to  Software  Requirements 
Specifications". 
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OVERVIEW  AND  COARSE  DECOMPOSITION 

OLTRIS  -  The  top  module  of  the  system 

1 .  Identification  OLTRIS 

2.  Type  Module. 

3.  Purpose  This  module  is  created  to  give  an  overview  over  the  total  system  and  its 

environment,  as  tniefly  outlined  in  Figure  1. 

The  system  is  part  of  a  research  project  with  the  aim  to  ”...  develop  a  data 
processing  concept  for  on-line  high  accuracy  point  positioning  ... 
systems".  To  find  the  best  solutions  for  this  concept,  there  are  pointed  out 
four  research  topics.  It  is  decided  that  the  concqrt  s^  be  materulized  in  an 
operational  system  for  use  in  aerial  photo-triangulation.  The  system 
should  woilt  on  digital  images. 


Figure  1;  The  On-Line  Triangulation  System  (X.TR1S  and  its  enyiroainent. 


4.  Function  The  system  will  perform  On-Line  Photo-triangulation.  Important  tasks  are; 

•  Initial  value  computation 

•  Sequential  and  simultaneous  bundle  adjustment 

•  Blunder  detection 

•  Self  calibration,  witii  test  ctf  significance  and  deteiminat^ty 
(•  Automatic  point  transfer) 

(*  Inclusion  of  Idnematic  GPS  data) 

Currently  the  program  "oltris"  is  working  as  an  experimental  program 
without  a  link  to  analytical  plotter  or  digital  workstation. 


5.  Subordinates  (See  Hgine  2)  This  module  is  composed  of  following  modules; 


USEIL.INT/XVIEW 

FLOW.OON 

-  User  interface  in  OpenWindows  (OpenLook) 

-  Admimstiation  module 

INPUT 

INTITAL 

ADJUST 

QUAI^OON 

REPORT 

-Data  input 

-  Initial  panuneter  values  compiitatiroi 

-  Bundle  adjustment 

-  Quali^  ccKUiol 

-  Ou^ut  of  results 

FUNC  UB 
MEMJVLLOC 

-  Library  for  cmnmon  functions 

-  Memory  allocation  functions 

CONTR  PAR 

MCXJEL  DATA 
E(2_SYST 

-  Control  parameters 

-  Data  for  entities  in  the  mathematical  model 

-  Normal  equatkm  system  and  related  data 

System 
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6.  Dqjendencies  A  session  with  diis  system  is  supposed  to  be  ctuiducted  interactively  by  a 
human  operator. 

Requiring  its  tespe^ve  source  equipmtmt,  data  input  is  siqiposed  to  be  any 
or  mote  of:  Digital  images,  cocmnnate  lists,  parameter  lists,  contrtd 
ccnnmands.  Input  may  occur  manually  or  from  disk,  depending  (m  the 
character  of  the  data. 

Ouqrut  of  results  is  sent  to  disk,  printer,  artd  display. 

8.  Resources  The  diri^  OD-line  plioto*triangiilation  system  is  to  be  implemented 
in  a  UraX  environnoent  on  a  Sun  workstation,  using  the  programming 
language  "C".  The  user  dialog  will  be  handled  throu^  a  wiiuiow  system 
(O^nVf^ndows). 

The  crqucity  (tf  die  system  should  be  practically  unlimited  concerning 
numbn  of  image  points,  object  points,  photognqrhs,  and  cameras. 

Response  times  should  te  well  within  the  limits  of  of  the  operatonr’s 
patience. 

Both  te^xmse  time  arid  memory  requirements  ate  to  be  more  clearly 
q>edfied. 

As  examples  of  possible  nemocy  lequiiemenis: 

•  For  normal  equadon  system  and  observation  data,  SO  photos  and  SOO  points,  each 
measured  in  4  photos,  may  need  2  Mb,  while  100  photos  and  1000  points  need  S  Mb. 
This  should  st^  resident  in  memory  during  seqneaiial  adjustment. 

•  For  digital  images  (or  puts  of  images)  available  for  on-screen  viewing,  10  image  patches 
of  S12xS12  8-bits  pixels  will  occiqty  2.5  Mb.  However,  only  a  few  need  be  in  memory  at 
onetime. 

•  For  digitized  p^ts,  which  are  defined  by  small  patches  ($  25x25  pixels).  1000  image 
points  may  require  riboot  1/2  Mb.  And  again,  not  all  are  needed  in  memory  at  one  time. 
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1 .  Introduction 

This  tWninv-nt  describes  the  decon^osition  of  the  On-Line  Triangulation  System  (OLTRIS), 
also  showing  dependencies  between  the  modules.  Hgure  1  gives  an  overview  over  die 
main  nwiriiiing  with  their  intenelatioas.  The  sht^  of  the  boxes  in  the  diagrams  denotes 
whedwr  the  modules  contain  mainly  functions  (rectangles)  or  primarily  ^ed  for  data  storage 
(toutKfod). 


Hgure  1;  An  overview  of  the  main  modnles  of  the  (XlltlS. 
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2.  USER_INT/XVIEW  -  User  interface  in  OpenWindows 

1.  Identification  USERJNT/XVIEW 

2.  Type  Module. 

3.  Purpose  Handle  all  dialogs  between  the  user  and  the  systeni,  which  theusercan 

lead  the  process  of  the  system.  Commands  are  to  be  given  menus, 

buttons,  and  ton  in  windows,  -  whatever  is  found  appropriate  for  the 
specific  case.  The  dialogs  include  also  any  manual  ii^t  of  parameters. 

4.  Function  In  a  OpenlA^ndowenviiooment  it  seems  natural  to  include  all  "set-up" 

activi^  for  dte  interactive  toc^  in  this  module,  while  action-taking  fimctions 
("Target  Functions")  belong  in  FLOW_(}ON. 

This  module  should  also  cootain  functions  needed  for  aunmunicating  witii 
Openl^findows  when  the  woridng  modules  need  to  change  or  add  elements 
of  the  user  interfiles. 

5.  Subordinates  This  module  contains  the  fitilowing  modules: 

ULDATA  -User  interface  for  all  data 

UI_S  YST  -  User  interface  for  the  system 

UI_CANV  -User  interface  for  the  canvases 

XV_NCyiHV  -  Notify  and  event  callbadc  functions 

6.  Dependencies  (See  Figure  1) 

•  Part  of  module  (X^TRIS 

•  Interacts  with  module 

FLOW_CX)N 

and  (possibly  indirectly)  widi 

INPUT 

INITIAL 

8.  Resources  Hardware:  Sun  workstations 

Software:  OpenMfindows  Version  3.0  ((3penLook) 

Libraries:  -  libxview.a 

-  lib(tigx.B 

-  libX.a 
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2.1.  U1_DATA  -  User  interface  for  the  data 

J.Idau^cation  ULPATA 

2.  Type  Module. 

3.  Purpose  Provide  panels,  windows  eic  for  tasks  of  data  handling. 

4.  Function  Conuius  the  lH  object  initieKntinnfiinctioMforihewindtwq^ic^tK*";  of 

data  iiqxit,  displaying  etc.  For  this  reason  the  objects  of  XView  data  types 
(popup  fireme.  menus,  buttons,  textpanels,  etc.)  win  be  created. 

5.  Subordinates  Ibis  module  contains  the  fi^owing  functions: 

CAMERA_UI  -iX^ndow  for  camen  data  input  and  di^lay 
IMACX_UI  •  Window  fior  image  data  input  and  display 
DISPLAY_UI  -V^ndow  for  for  di^layttf  data  and  results 
INOBS.UI  -  >K^ndow  for  input  at  observations 

6.  Dependencies  'Pan  of  module  USEIL.INT/XVIEW 

*  Interacts  widi  die  nwhiif- 

FLOW.CON 

8.  Resources  Software:  OpenWindows  Version  3.0  (OpenLocA) 
libruies:  •  libxview.a,  libolgx.a.  UbX.a 

2.2.  UI_SYST  -  User  interface  for  the  system 

1 .  Identification  UI_SYST 

2.  Type  Module. 

3.  Purpose  Provides  windows  for  processing  of  (X.TRIS. 

4.  Function  Contains  the  UI  object  initialiTaiion  functions  for  the  window  qqdications  of 

die  system  processing  etc.  For  diis  reason  the  otgects  of  XView  data  types 
(base  and  popqi  canvase,  menus,  buttons,  textpanels,  etc.)  will  be 

created. 

5.  Subordinates  This  module  contains  the  following  foncdons  (windows): 

OLHPOLUI  -Dialog  for  program  flow  oontiol 

(XTMAIN.UI  -Dialog  for  adgukmentoontitd  parameter 

PARASET.UI  -Dialog  for  initial  value  control  parameter 

CENSTAT.UI  -  Kalog  for  program  flow  control 
l*D>iFO_UI  -Dialog  for  program  flow  control 

(2UALOON_UI  -Diahy  for  athament  control  parameter 

REPCXIT.UI  -  Dialog  for  mitialvtthie  control  parameter 

6.  Dependencies  •Pattofmodule  U5ERJNT 

*  htteracta  widi  the  module 

FLOW.CXW 

8.  Resources  Software:  OpenWndows  Veraon  3.0  (OpenLoolO 
Libraries:  -  Hboiview Ji,  libtdgxj,  UbXji 
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2.3.  UI_CANV  •  User  interface  for  the  canvas  applications 

l.Idaa^icatioH  ULCANV 

2.  Type  Module. 

3.  Purpose  Provide  canvases  for  OLTKIS  .or  diq>layine  positioas  of  measured  image 

pointt  and  images  in  the  pborogiaimDeiric  bac^ 

For  image  coerdinates:  In  the  final  s^tem,  there  should  be  one  or  more 
canvases  for  diqtlaying  (parts  of)  digital  images.  Positioning  on  the  images 
is  to  be  done  by  mnuMt  cursor.  (Image  matching  techni^ies  will  also  be 
used  for  precise  positioning,  but  therefore  an  inteifiice  must  be  provided  to 
liiA  tte  oh  system  with  a  tenylatc  matching  program.) 

4.  Function  Contains  the  UI  object  inidalizatioo  functions  for  die  canvas  qiplications  of 

the  system.  For  this  reason  the  objects  of  XView  data  types  (canvase, 
menus,  buttons,  etc.)  wiU  be  created. 

5.  Subordinaus  This  module  contains  the  following  functions  (canvases): 

CANV _BLOCX_UI  -  Canvas  for  die  photogiammetric  block. 

CANVJDdG.UI  -  Canvas  for  diqilay  of  images  and  meas.  points 

6.  Dependencies  •Part  of  module  USER^Ihn' 

•  Interacts  widi  die  mothde 

INPUT 

8.  Resources  Software:  OpenVfindows  Version  3.0  (OpenLook) 

Libniy:  •  libsview.a.  libolgx.a.  libX.a 

2.4.  XV_NOTEV  •  Notify  and  event  callback  functions 

1.  Identification  XVJfOTEV 

2.  Type  Module. 

3.  Purpose  Provide  notify  and  event  callback  functkais. 

4.  Function  Window  system:  Set  iq>  the  user  dialog  for  object  parameter  iiqMiL 

5.  Subordinates  This  module  is  composed  of  the  functions: 

DLJMACS  -Dialog  for  image  data  irqmt 
IX^CAMERA  -Dialog  for  camera  data  input 
DL,_MCAM  -Dialog  for  manual  camera  daa  input 

IX^OPC  'Dialog  for  object  point  cooidinaies  input 

DL_FOP  -Dialog  for  exterior  ofieauuion  parameters  input 

6.  Dependencies  •  Part  of  module  USER_INT 

•  Interacts  with  the  module 

INPUT 

8.  Resources  Software:  OpenWindows  Versioa  3.0  (OpenLook) 

Ufanries:  -  tibxviewa,  ]ibolgx.a.  UbXa 


3.  FLOW_CON  •  Administration  module 

1.  Idadfiattion  FLOW_CON 

2.  Type  Module. 

3.  Purpose  To  peifixm  the  cential  control  of  the  system,  based  on  user  commands  and 

intermediate  results  in  the  system. 

4.  Function  This  module  contains  the  main  piognun  and  functions  ti^iich,  based  on  user 

commands  and  intermediate  re^ts,  will  administer  the  contrd  flow.  All 
calls  to  tnain  functions  in  the  ’worl^g”  modules  of  the  system  should  then 
come  ficm  this  module. 

When  using  a  window  system,  this  module  includes  the  Target  Functions" 
activated  ^  interacdve  tods  (OpenWindows)  to  interpret  the  events  and 
decide  actKMts  to  be  taken. 

5.  Subordinates  This  module  contains  the  fdlowing  modules: 

FC.MAIN  -  Tcp-level  flow  contid 

FC_MEASURE  -  Interfisce  to  the  INPUT  module  (Measurement) 

FC.INPUT  -  Interface  to  the  INPUT  module 

FC.INCLUDE  -  Inclusion  at  observations 

FC_EXCLUI£  -  Exclusion  of  observations 

PC_PISPLAY  •  Di^lay  of  datt  and  results  on  screen 

FC_UPDATE  •  Interfax  to  the  ADJUST  module 

FC.QUALC^  -Interne  to  the  (^AL_C(^  module 
FC_WaPORT  -  Inter&ce  to  the  REPCHIT  module 

6.  Dqiendencies  (SeeHgurel) 

‘IspartrfOLTRB 

*  Supposed  to  interact  widi  all  the  main  modules  conqxrsing  OLTRIS 
8.  Resources  None  required. 
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3.1.  FC_MAIN  -  Top-level  flow  control 

/.  Idendficcuion  FC_MAIN 

2.  Type  Module. 

3.  Purpose  Provide  the  central  admimstration  of  OLTRIS. 

4.  Function  This  module  performs  the  initializadoa  of  X^ew,  of  the  user  interface  and 

OLTRIS  coaqxments 

5.  Subordinates  This  module  contains  the  function: 

XC%.T  -  Top-level  flow  control 

6.  Dependencies  •  Is  part  module  FLOW_OON 

*  Interacts  with  the  modules 
USER_INT/XVIEW 

PC_>ffiASURE 

FC_INPUT 

FC.INCLUDE 

FC^CLUDE 

FC.DISPLAY 

FC_ADJUST 

FC.QUALCON 

FC_REP(»T 

8.  Resources  None. 


3.2.  FC_MEASURE  •  Interface  to  the  INPUT  module 

1.  Identification  FC_MEASURE 

2.  Type  Module. 

3.  Purpose  Provides  an  intei&ce  to  dw  nxxlule  INPUT  (Measurement). 

4.  Function  Hus  module  peifbnns  detailed  cx»trol  of  the  opei^ukuis  in  the  INPUT 

module  with  req>ect  to  input  of  image  point  condinates. 

5.  Subordinates  This  module  contains  the  function: 

NEXT_01PC  -  Iiqttit  of  next  image  point  observaticms 

M0RE_01PC  -  Input  of  mme  than  one  image  point  observaticxis 

SI^C.OIPC  -  Input  of  a  qiecific  image  point  observations 

MAN.OIPC  -  Ntoual  input  of  image  pcwt  observaticxis 
TIE_OIPC  -  Input  of  image  point  oteeivations  of  tie  points 

CONPO.OIPC  -  Input  of  image  point  observations  of  control  pc^ts 

6. Dependeru:ies  •  Is  part  of  module  FLOW_CON 

•  Interacts  widi  die  modules 

UI_SYST 

INPUT 

8.  Resources  None  requiiecL 

3.3.  FC_INPUT  -  Interface  to  the  INPUT  module 

1.  Identification  FC_INPUT 

2.  Type  Module. 

3.  Purpose  Provides  an  interface  to  the  module  INPUT. 

4.  Function  This  mobile  performs  detailed  control  of  die  operations  in  die  INPUT 

mcxhile. 

5.  Subordinates  This  module  contains  the  following  fuiKXkxis: 

XV.CAMERA  -  Input  d  camera  calibratkxi  data 
XV_IMAGE  •  Ii^  image  data 

XVjOOPC  -Injiut  of  object  point  (diservations 

XVJOEOP  -Input  of  exterior  orientation  parameters 

6.  Dependencies  •  Is  part  of  module  FLOW.CON 

•  Interacts  with  die  modules 

ULPATA 

UI_SYST 

INPUT 

8.  Resources  None  required. 


3.4.  FC_INCLUDE  -  Inclusion  of  observations 

1 .  Ident^cation  FC_INLUDE 

2.  Type  Module. 

3.  Purpose  Includes  (additional)  observatims  in  the  adjustmenL 

4.  Function  This  module  includes  additx»al  observations  in  the  adjustment 

5.  Subordinates  This  module  contains  the  functions; 

ALL_IN(Xi_OIPC  -  Include  of  all  unused  ipc  observations 

ALI._INCL_(XX*C  -  Include  ctf  all  unsused  opc  observations 
AIJL_INCL_OEOP  -  Include  of  all  unsused  eop  observaticms 

SPEC_lNCL_OIPC  -  Include  of  one  qtecific  unused  ipc  observations 

SPEC_INCL_OOPC  -  Include  of  one  q)ecific  unused  opc  obs. 

SI%C_IN^_OEOP  -  Include  of  one  specific  unused  e^  obs. 

6.  Dependencies  •  Is  part  module  FLOW_CON 

*  Interacts  widi  the  modules 

INPUT 

8.  Resources  None  required. 

3.5.  FC_EXCLUDE  •  Exclusion  of  observations 

1 .  Ident^ication  FCJEXCLUDE 

2.  Type  Module. 

3.  Purpose  Excludes  obsorvations  from  the  adjustmenL 

4.  Function  This  module  excludes  observations  from  dbe  adjustment 

5.  Subordinates  This  module  contains  the  functions: 

SPEC_EXCL_01PC  -  Exclude  of  one  specific  unused  qx;  obs. 

SPEC_EXCL_OOPC  -  Exclude  of  one  specific  unused  opc  obs. 
SPEC_EXCL_OE(X’  -  Exclude  of  me  specific  unused  e^  obs. 

6.  Dependencies  *  Is  part  of  module  FLOW.CON 

*  biteracts  with  the  modules 

<3UAL_CON 

8.  Resources  None  requited. 


3.6.  FC_DISPLAY  •  Display  of  data  and  results  on  screen 

1 .  Idmt^cation  FC_DISPLAY 

2.  Type  Module. 

3.  Purpose  Rt>vide  an  intetface  to  the  module  REPORT. 

4.  Function  This  module  perfonns  detailed  contnd  of  the  operations  in  displaying  data 

and  results  oa  screen. 

5.  Subordinates  This  module  is  coiiq>osed  d  the  modules: 

XV.OIPC  -  Display  of  ^  observations  in  measured  sequence 

XV_OIPC_IMG  -  Di^lay  of  ipc  observadtms  m  image  sequence 
XV_01PC_(X^  -  Di^lay  of  ipc  observatitMis  in  objea  point  sequence 

XV_OOPC  -  Display  of  object  point  coordinate  obxrvaticHis 

XWJOEGP  •Diq>lay  of  exterior  orientation  parameter  observ. 
XV_IMAGES  -IMsplay  of  image  data 

XVjyPT  -  Di^Iay  d  objM  point  cordinates 

XV_CAMERA  -  Di^lay  of  camera  calibration  data 

XV_SOLEOP  -  Di^lay  (rf  cop  results 

XV_SOLOTC  -  Display  d  opc  results 

XV_SOLRESID  -  Diqday  of  residuals 

XV_SOLN(%MAL-  Di^lay  of  normal  equation  system 

6.  Dependencies  •  Is  part  of  module  FLOW.CON 

•  Interacts  with  the  modules 

ADJUST 

REPC»T 

8.  Resources  None  required. 

3.7.  FC_UPDATE  ■  Interface  to  the  ADJUST  module 

1 .  Identification  FC.UPDATE 

2.  Type  Module. 

3.  Purpose  Provides  an  interface  to  the  module  ADJUST. 

4.  Function  This  module  performs  control  of  the  (iterations  in  the  ADJUST  module 

(qtecially  in  the  function  ADJtELlN). 

5.  Subordinates  This  module  is  composed  d  the  functi(»: 

FC_SEQ  -  Update  secpientially 

FC.SIM  -  Update  and  relineaiize  simultaneously 

6.  Dependencies  •Is  part  of  module  FLOW.CON 

•  Interacts  with  the  module 

ADJUST 

8.  Resources  None  required 
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3.8.  FC_QUALCON  -  Interface  to  the  QUAL_CON  module 

1.  Identification  FC_QUALCON 

2.  Type  Module. 

3.  Purpose  Provides  an  inteiface  to  die  module  QUAL_CON. 

4.  Function  This  module  performs  detailed  control  of  the  operations  in  the  QUAL_CON 

module. 

5.  Subordinates  This  module  contains  the  following  fimctions: 

BLUN_DET  -  Blunder  detection 

SELF_CALIB  -  Self  calibraion 

6.  Dependencies  •  Is  part  of  module  FLOW_CON 

•  Interacts  widi  the  modules 

QUAL_CON 

8.  Resources  None  required. 

3.9.  FC_REPORT  •  Interface  to  the  REPORT  module 

1 .  Identification  FC_REPORT 

2.  Type  Module. 

3.  Purpose  Provides  an  interface  to  the  module  REPORT. 

4.  Function  This  module  performs  detailed  control  erf  the  operations  in  the  REPORT 

module. 

5.  Subordinates  This  module  contains  the  following  functions; 

REP_DISK  -  Flow  control  for  the  saving  data  on  disk 
REP_PRINT  -  Flow  control  for  the  printing  of  reports. 

REP_PLX)T  -  Flow  control  for  the  plotting  of  dara 

6.  Dependencies  •  Is  part  of  module  FLOW_CX)N 

•  Interacts  with  the  module 

REPCSIT 

8.  Resources  None  required. 


4.  INPUT  -  Data  input 

1.  Identification  INPUT 

2.  Type  .  Module. 

3.  Purpose  Receives  iiq>ut  of  data  firom  outside  the  system  and  from  measurements  on 

digital  images  (Digital  Photogiammetric  Woriatation),  and  organizes  and 
stores  it  in  internal  data  structures. 

4.  Function  Funcdoas  of  this  module  will  control  the  input  of  data,  and  organ^  and 

store  it  in  internal  data  structures.  Data  may  come  firom  the  user  dialog,  fiom 
files,  or.dir«^y  as  measurements. 

5.  Subordinates  This  module  is  cmnposed  ctf  the  modules  (See  figure  2): 

IN_IMAGE  •  Input  of  digital  images 

IN_CON_P  AR  -  Input  of  control  parameters 

IN_OBJ_PAR  -Input  of  parameters  for  objects 

IN_OBSERV  -Ir^ut  of  observation  data 

6.  Dependencies  (See  Hgure  1) 

•  The  rtxxlule  is  part  of  nxxiule  OLTRIS 

•  The  rtxxlule  interacts  with  the  modules 

FLOW_CON 

USER.INT 

DIG_IM 

CONTR_PAR 

MOC®LJ5ATA 

8.  Resources  Hardware :  Disk 

Software :  Interface  to  external  sources. 


Figve  2:  Overview  of  the  INPUT  nxxiule. 
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4.1.  IN_IMAGE  •  Input  of  data  from  digital  inures 

1 .  Identification  IN.IMAOE 

2.  Type  Module. 

3.  Purpose  Handles  the  input  of  data  from  digital  images. 

4.  Function  This  module  contains  fimctioosnecessaiy  to  read  and  store  digital  imi^ 

data  (image  coonUnate  observatims)  in  the  system's  internal  storage. 

(The  first  version  of  the  input  of  digital  images  data  into  the  Systran  will  use 
data  sets  consisting  of  lists  of  image  coonlinates  measured  on  adigiral 
wodcstation  An  extraided  version  should  provide  an  intraface  to  an  «miplati» 
matching  program  to  organize  die  intermediate  storage  of  this  data  in 
common  data  structures.) 

5.  Subordinates  This  module  is  composed  d  the  modules  (functions): 

...  (Not  yet  considered  and  installed) 

6.  Dependencies  (SeeHgure2) 

•  Is  part  module  INPUT 

•  Interacts  widi  die  modules 

FLOW.OON 

USER_INT/XVIEW 

8.  Resources  Software:  biteiface  to  a  template  matching  program. 

Hardware:  Digiral  photogrammetiic  wortetation 

4.2.  IN_CON_PAR  -  Input  of  control  parameters 

1 .  Identification  IN_CON_PAR 

2.  Type  Module. 

3.  Purpose  Handles  the  input  d  control  parameters. 

4.  Function  This  module  contains  funcdons  fix’  fetching,  organizing,  and  storing  the 

control  parameters.  Data  may  come  fiom  ™niii»l  input  or  fiom  file. 

: .  Subordinates  This  module  contains  the  following  functions: 

...  (Not  yet  considered,  because  the  input  of  craittol  parameter  comes 
direcdy  mm  file  or  is  set  in  Module  USER_1NT/XV1EW) 

6.  Dependencies  (SeeHgure2) 

•Is  part  of  module  INPUT 

•  hiteracts  with  the  modules 

FLOW_CX)N 

USERJNT/XVIEW 

CX)NTRJ>AR 

8.  Resources  None  required. 


4.3.  IN_OBJ_PAR  -  Input  of  parameters  for  objects 

l.Ident^cation  IN_OBJ_PAR 

2.  Type  Moduk. 

3.  Purpose  Handle  the  input  of  parameters  for  otgects. 

4.  Function  This  module  ccmtains  functions  for  fetching,  organizing,  and  storing 

parameters  for  various  objects  in  the  data  stnictuies  of  the  system.  Data  may 
come  fiom  manual  input  or  fiom  file. 

5.  Subordinates  This  module  is  composed  of  the  following.functions: 

IN.OEOP  -  Input  exterior  orientation  parameters 

Df.OOPC  -  bput  otgect  point  coordinates 

IN_CAMERA  -  li^ut  camera  parameters 
IN_IMG  -  Ii^ut  image  dm 

6.  Dependencies  (SeeHgure2) 

•  Is  part  module  INPUT 

•  Interacts  with  the  modules 

FLOW.CON 

USER_INT 

M(X)EL_DATA 

8.  Resources  Not  yet  considered. 

4.4.  IN_OBSERV  -  Input  of  observation  data 

1 .  Identification  IN_OBSERV 

2.  Type  Module. 

3.  Purpose  Handle  the  input  of  observatioo  data. 

4.  Function  This  module  contains  functions  for  fetching,  organizing,  and  storing 

observation  data  in  the  data  structures  ^  the  system.  may  crane  fiom 

manual  itqnit  or  fiom  file. 

5.  Subordinates  This  module  is  conq)osed  of  the  fiuxxions: 

IN.OIPC  -  Iiqrut  im^  point  cooidinates  observations 

IN.OOPC  -  Input  djject  poim  cooidinates  observations 

IN.QEOP  -  hiinit  exterior  orientation  parameters  observations 

6.  Dependencies  (See  Figure  2) 

•  Is  part  of  module  INPUT 

•  Interacts  with  die  modules 

FLOW.OON 

USERJNT 

MODELJ^ATA 

8.  Resources  Not  yet  consideted. 
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5.  INITIAL  •  Initial  parameter  values  computation 

1.  Identification  INITIAL 

2.  Type  Moduk. 

3.  Purpose  Conqnite  dte  initial  values  for  aU  unknown  paranieters.  Necessary  because 

of  the  non-linear  equation  system. 

4.  Function  Rovide  initial  values  of  exterior  orientation  parameters  aiKl(^  object  point 

coordinates  for  an  aerial  Mock  using  diffetent  methods  of  conqtuting: 
intersection,  resection  in  qtace  and  relative  orientation. 

5.  Subordinates  This  module  is  cooqxMed  of  the  modules: 

INTT.ECX’  -Initial  values  for  exterior  orientation  parameters 

INIT.OPC  -hiitial  values  for  otgect  point  coordinates 

6.  Dependencies  •  The  module  is  part  module  (X,TRIS 

•  The  module  interacts  with  the  modules 

FLOW_CX)N 

FUNC_LIB 

CONTK-PAR 

MOC«L_DATA 

8.  Resources  None  required. 

5.1.  INrr_EOP  •  Initial  values  for  exterior  orientation 

1.  Identyication  INITJBOP 

2.  Type  Module. 

3.  Purpose  Provide  initial  values  for  exterior  orientation  parameters. 

4.  Function  This  module  contains  functions  for  checking  the  presence  of,  and  if 

necessary,  ccxiqtuting  (or  otherwise  provide)  initial  values  for  exteriOT 
orientation  parameters. 

5.  Subordinates  This  module  is  ooriqiosed  of  the  two  modules: 

RELQR  Functions  for  idative  orientation 

RESEC  Functions  for  resection  in  qxKe 

6.  Dependencies  <  Is  part  of  module  INITIAL 

•  Interacts  with  die  modules 

FLOW.OON 

FUNC_LIB 

OCWTRJPAR 

MODEL_PATA 

8.  Resources  None  requited. 


Docjto. 

rate 

SD.03 

ition  Descripnem 

Ver.  Page  17 
2.1  of  41 


5.1.1.  RELOR  -  Relative  orientation 

1 .  Identfflcation  RELXXt 

2.  Type  Kfodiik. 

3.  Purpose  CanqwdBilie  initial  values  of  exterior  crientatknpaiainelas  by  lebttive 

nrit-nwtinn  by  using  die  following  method:  coplmari^  equations,  sin^e 
modeL 

4.  Function  Provide  initial  values  of  exterior  orientation  paiameters  and  put  it  into  the 

data  stiuctutes 

5.  Subordinates  This  module  is  cooqiosed  of  the  functions: 

INTT  JtELCK  •  checks  some  conditions  for  relative  (mentation 

M(X)_PEL(Xt  -relative  orientation  for  models 

LOCEOP_STAR.TVAL  -  provides  start  values  to  conqxite  relor 

6.  Dependencies  *  The  module  is  part  c£  module  INTr_EOP 

•  The  module  interacts  widi  the  modules 

FLOWjCXlN 

FUNC.LIB 

CC»1TR^PAR 

M(XeLJ3ATA 

8.  Resources  None  required. 

5.1.2.  RESEC  •  Resection  in  space 

1 .  Identffication  RESEC 

2.  Type  Module. 

3.  Purpose  lYovide  initial  values  for  exterior  orientation  parameters. 

4.  Function  This  module  contains  functions  for  checking  ocxiditions  of  the  resection, 

ptqwting  data  for  the  resection  and  executing  the  resection  in  qtace  by  an 
adjustment 

5.  Subordinates  This  mochile  is  oooqxtsed  of  the  following  functions: 

lNrr_ECX*  checks  the  conditions  for  resection  in  yce 

INIT_ABSE(X*  contains  some  functions  for  configuration  tests 

INrr_RESEC  data  prqwraiion  for  resection  in  qtace 
RES_AF1t.OX^VAL  provides  start  values  for  resectiem  in  qiaoe 
RES_COMP_RESK  atyustment  of  resection  in  space 

RES_PARJK  conpuies  the  partial  derivatives  of  observations 

6.  Dependencies  •  b  part  of  module  INrr_E(^ 

•  Inteticts  with  the  modules 

FUNC_IJB 

a»mu»AR 

MGOELJDATA 

8.  Resources  None  requited 


6.  ADJUST  -  Bundle  adUustment 

l.Ident^ication  ADJUST 

2.  Type  Module. 

3.  Purpose  Pofbnn  die  bundle  adjustment  tasks,  shnultaneous  as  well  as  sequendaL 

4.  Function  Build  and  fiKtorize  nonnal  equations,  iqidate  fiKtorized  noimal  equations 

sequential^,  compute  solution  vector,  and  update  die  parameter  vector. 

5.  Subordinates  This  module  is  composed  oftbefidlowing  modules  (See  figure  3): 

BUILD  -  Build  nomud  equations 

CHFAC  -  Factorize  normal  equations 

lUDAlB  -  Update  Smtorized  normal  equation  system 

SOL^  -  Solve  normal  equations 

PARAM  -Update  object  parameters  (unknowns) 

6.  Dependencies  (See  Hgute  1) 

•  The  module  is  part  of  module  (X.TRIS 

•  The  module  interacts  with  the  modules 

FLOW_CON 

FUNC_LIB 

C(»mL.PAR 

EQ_SYST 

MODELJ^ATA 

8.  Resources  None  required. 


PlfM  3:  Owriew  of  dw  ADIUST  andala. 


Docjto.  Title 
SD.03  Deg 


ioo 


yer.  Page  20 
2.1  of  41 


System  I 
OLTRISl 


ompQgit 


Deacription 


Authais) 

Dme 

Ho.tk 

15.03.1992 

6.1.  BUILD  •  BuUd  normal  equations 

l.Identificatum  BUILD 

2.  Type  Module. 

3.  Purpose  Build  the  normal  equations  (A'H'A,  A'n*l). 

4.  Function  Go  through  aD  observadoo  dtta  and  coaq>ute  and  add  their  conoibutions  to 

the  normal  equadoo  system. 

5.  Subordinates  This  module  is  oonqtosed  of  the  fiidlowing  functions; 

Add  contributions  to  tile  normals 
SIM-BUILD  -  One  coDoplete  observation  equation 

SIM_BLD_IPC  -  Observations  of  one  image  point 

SIMJBLO.OPC  -  Observations  of  obj.  cot^  of  erne  obj.  point 
SIM_BLD_E(X*  -  Observations  ctf  ext  ori.  param.  of  one  image 

6.  Dependencies  (SeeRgure3) 

•  The  modute  is  part  cf  module  ADJUST 

•  The  module  interacts  with  the  modules 

FLOW_CX)N 

FUNC_LIB 

CONTR^PAR 

EQ_SYST 

MCK5ELJ3ATA 

8.  Resources  None  required. 

6.2.  CHFAC  •  Factorize  normal  equations 

1 .  Ident^ication  CHFAC 

2.  Type  Module. 

3.  Purpose  Factorize  the  normal  equation  system,  Le.  form  an  upper  triangular  matrix. 

4.  Function  The  normal  equation  system  is  foctorized  by  Cbolesky  fimtorization. 

5.  Subordinates  This  module  is  oonyosed  of  the  following  functions: 

SIM_CHFAC  -  Oxtiesky  foctorization  (general) 

SIM_CHFAC_HYP  •Oxti.facL  (qiecial  hyper-row  treatment) 

SIM_CHFAC_RES  -Cheti.  foct  (resident  pan  of  the  normals) 

6.  Dependencies  (See  Figure  3) 

•  The  module  is  pan  of  module  ADJUST 

•  The  module  interacts  with  the  modules 

FLOWJXW 

(CONTILPAR) 

EQ_SYST 

8.  Resources  None  required. 
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6.3.  UPDATE  -  Update  factorized  normal  equation  system 

1 .  Ident^ication  IIFDAIE 

2.  Type  Module. 

3.  Purpose  Uixbae  die  fiKtorized  mxmal  equation  system  wiA  the  cootribotion  &om  one 

or  moR  additional  observations. 

If  "continuous  le-Hneaiuation"  is  inqdanented,  a  module  should  be  added, 
sdikh  will  organiie  die  "simultaneous"  down-  and  iqidate  widi  old  and  new 
contributions  from  a  qiecified  set  of  observations. 

4.  Function  F4tch  the  qiecified  obaervation(s),  compute  the  observation  equation 

coefikien^  and  perform  the  upi^. 

An  algorithm  based  rai  Givens  lYansCDmiation  is  used  for  the  iqxlate. 

5.  SubordUnates  This  module  is  conqposed  of  the  functions: 

UDJPC  -  Update  one  ima^  point  observation 

UD_1IPC  -  Organize  ipdate  with  cme  image  point  observatimi 

UD_SQIPC  -  Organize  iqidaies  with  a  sequence  of  image  points 

UD_CMPC  -O^anizeiqidates  with  image  points  concerning  one 

obgKt  point 

UD_IMIPC  -O^uuze  updates  with  sequence  of  image  point 
observatkMis  from  one  image 

UD.OPC  -  Update  by  observations  frm  coe  object  point 

UD.IQPC  -Or^unzeipdate  by  observations  from  one  Direct  point 

UD_EOP  -Ui^tate  by  observiuions  of  exterior  orientation 

UD.IEOP  •  Organize  iqxlate  by  observations  of  exterior  orientation 

GT.UFDATB  -  Perform  Givens  TYansformation  iqidate 
GT_PODATB  -  Perform  Givens  Transformation  downdate 

6.  Dependencies  (See  Figure  3) 

•  Is  part  module  ADJUST 

•  Interacts  widi  the  module 

FLOWJDON 

FUNC_LIB 

CONTIU*AR 

E(i.SYST 

MODELJ5ATA 

8.  Resources  None  required. 
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6.4.  SOLVE  -  Solve  normal  equations 

1.  Idendfication  S(H.VE 

2.  Type  Moduk. 

3.  Purpose  S<dve  die  normal  equatioa  system,  i.e.  compute  die  conecdons  to  the 

cuncndy  valid  values  of  die  parameters. 

It  should  be  possible  to  solve  for  a  specified  subset  of  die  parametm, 
namely  exterior  oientatkm  paxameteis  fir  die  last  image(s),  and  ot^ 
comdiiuues  for  specified  points. 

4.  Function  The  mxxnal  equation  system  is  solved  by  back  substitution. 

5.  Subordinates  This  module  is  conqiosed  of  the  fcdlowing  modules: 

SOL_BACKSUB  -  Back  substitution 

S(X._E(^  -  Back  substitution  for  ext  on.  parameteis 

SOL.OPC  -Back  substitution  for  object  point  comdinates 

6.  Dependencies  (See  Figure  3) 

•The  module  is  pan  of  module  ADJUST 

•  The  module  interacts  with  the  modules 

FLOW_CON 

(CONTRJ>AR) 

ECLSYST 

8.  Resources  N<»e  required. 

6.5.  PAR  AM  •  Update  object  parameters  (unknowns) 

1 .  Idendfication  PARAM 

2.  Type  Module. 

3.  Purpose  Update  the  values  of  die  unknown  parameters. 

4.  Function  The  solution  vector  of  the  equation  system  is  added  to  the  currmt  parameter 

vector.  It  should  be  checked  that  the  solution  vector  is  computed  after  the 
most  recent  iqxlate  of  the  equation  system. 

Later  versions  of  this  module  should  be  aUe  to  update  only  a  qiecified 
subset  of  the  parameters.  Normally  that  will  apply  to  parameters  for  recendy 
included  ima^s)  and  object  points,  which  are  likely  to  have  large 
conecticms. 

5.  Subordinates  This  module  is  oon^osed  of  the  following  ftmctions: 

AD_RELIN  -  Simultaneous  adjustment 

AD_PAR^UFDATE  -Update  of  the  parameter  vector  (unknowns) 

6.  Dependencies  (See  Figure  3) 

•  The  module  is  pan  of  module  ADJUST 

•  The  module  interacts  with  the  motfades 

FLOW_(X)N 

(C»NTRJ»AR) 

EQ.SYST 

MOC«L_DATA 

8.  Resources  None  required. 
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7.  QUAL_CON  -  Quality  control 

Uden^icadon  QUAL_GON 

2.  Type  Module. 

3.  Purpose  Evaluates  the  coupled  results,  the  observations,  and  the  mathematical 

model  used. 

4.  Function  Conqnites  some  general  statistics,  characterizing  the  results  as  a  whole, 

perfcms  Uunder  detectum,  and  checks  dctenninability  and  significance  of 
additional  parameters. 

5.  Subordinates  This  module  is  ctKi^xised  of  the  modules: 

GEN_STAT  -  General  statistics 

BLUNDER  -Blunder  detection 

AP_SIGNIF  -  Signifimnce  of  additional  parameters 

AP_DETERM  -  Determinability  of  additional  parameters 

6.  Dependencies  (See  Figure  1) 

•  The  module  is  part  of  module  (X.TRIS 

•  The  moduk  interacts  with  the  modules 

FLOW^CON 

CONTR_PAR 

E(i.SYST 

MOMLJ3ATA 

8.  Resources  Not  yet  considered. 

7.1.  GEN_STAT  -  General  statistics 

1.  Ident^cation  GEN_STAT 

2.  Type  Module 

3.  Purpose  Confutes  general  statistics  from  the  adjustment 

4.  Function  Conqiute  VTPV,  op,  redundancy,  etc. 

5.  Subordinates  This  module  contains  the  following  function: 

(yi_GENSTAT  -  CoDqwtation  d  general  statistics 

6.  Dependencies  •The  module  is  part  of  module  QUAL_CX>N 

•  The  mothile  interacts  with  the  modules 
FLOW_C»N 

EQ_SYST 
ADJUST 

None  required. 
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7.2.  BLUNDER  -  Blunder  detection 

1.  Identification  BLUNDER 

2.  Type  Module 

3.  Purpose  Detects,  localizes  and  initiates  removal  of  possible  blunders. 

4.  Function  Test  all  observations  for  the  presence  of  a  gross  error  widi  the  meduxl  of 

Batuxla's  dtda-snooping.  versimi  is  including .  Later,  an  extended  version 
should  testa  selec&«ofolMervationstosavecoiiq>atatk»time. 

Along  with  research:  Possible  improvements  or  variants  of  data-snocqring, 
as  well  as  odKT  medKxls  will  be  inq>lemented  and  tested. 

5.  Subordinates  This  module  contains  the  fcdlowing  functions: 

QC.SNOOP  -  Data  snooping  of  all  observations 

QC.EXCLUDEBLUNDER  -  Exclude  of  one  obs.  from  the  nonnals 

6.  Dependencies  •  The  module  is  part  of  module  QUAL_CON 

•  The  module  interacts  with  the  modules 

FLOW.CON 

CONTR_PAR 

EQ_SYST 

MC»EL,_DATA. 

8.  Resources  None  required. 

7.3.  AP_S1GNIF  •  SigniHcance  of  additional  parameters 

1.  Identification  AP.SIGNIF 

2.  Type  Module 

3.  Purpose  Test  the  significance  of  the  computed  additional  parameteis,  and  initiate 

recalctdatioa  without  die  exclutkd  parameters. 

4.  Function  Not  yet  described. 

5.  Subordinates  Not  yet  deconqxKcdAlso  this  module  may  be  con^xxed  of  alternative 

modules,  during  research  and  test  of  methods. 

6.  Dependencies  •  The  module  is  pan  of  module  QUAL_CX}N 

•  The  module  intenurts  with  die  module 

FLOW.CON 

cx»rrR..PAR 

EQ.SYST 

MODEL_DATA 

8.  Resources  Not  yet  considered. 
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7.4.  AP_DETERM  •  Determinability  of  additional  parameters 

1 .  Idatt^catioH  AP_Pt'ilikM 

2.  Type  Module 

3.  Purpose  Test  the  detenmnalnli^  of  the  additranalpanmeters,  and  initiate 

recalculation  without  im  excluded  parameten. 

4.  Function  Not  yet  described. 

5.  Subordinates  Not  yet  decooqxjsedAlso  this  module  may  be  conqxised  of  alternative 

modules,  during  research  and  test  of  metlnds. 

6.  Dependencies  •  The  module  is  part  of  module  QUAL_CON 

•  The  module  interacts  with  the  modules 

FLOW.CON 

CXWTILPAR 

EQ_SYST 

MOraUJATA 

8.  Resources  Not  yet  omsidered. 

8.  REPORT  •  Output  of  results 

1.  Identyication  REPORT 

2.  Type  Module. 

3.  Purpose  Format  and  ouq>ut  results;  for  displaying  on  the  screen,  as  repons  on  paper, 

on  secondary  storage  for  further  use.  ITotting  is  included,  on  display  (and 
on  paper). 

4.  Function  Separate  modules  will  oigai^  the  handling  of  results.  Display  it  on  screen, 

print  it  on  pq>er  and  store  it  in  files  on  disk. 

5.  Subordinates  This  module  is  cooposed  of  the  fcdlowing  modules; 

RE_D1SF1AY  -  Di^lay  results  on  the  screen 
RE_PRINT  -  noduce  report  on  paper 
RE_STORE  -  Store  results  on  disk 

6.  Dependencies  (See  Rgure  1) 

•  The  module  is  part  of  module  OLTRIS 

«  The  module  interacts  with  the  modules 

FIjOW_CX>N 

CX»frR_PAR 
MOi:£LJ>ATA 

Hardware:  -laser printer 

-  matrix  printer 
-disk 

-  screen  of  a  workstation  or  tetminal 


8.  Resources 


8.1  RE_DISPLAY  -  Display  of  data  and  results  on  screen 

1 .  Identification  RE_DISPLAY 

2.  Type  Module. 

3.  Purpose  Display  of  input  data  (observations  and  objects),  of  intermediate  results  and 

the  solution  on  the  screen . 

4.  Function  The  displaying  of  the  difforent  data  will  be  excuted  in  different  functions . 

5.  Subordinates  This  module  is  conqmsed  of  the  following  modules: 

XV_DSP_01PC  -  Display  of  ipc  obs.  in  measured  sequence 

XV_DSP_OIPC_IMG  -  Display  of  ipc  observations  in  image  sequence 

XV_DSP_OIPC_OPT  -  Di^lay  of  ipc  obs.  in  object  point  sequence 

XV_DSP_OCH*C  -  Display  of  object  point  coond^ate  observations 

XV_DSP_OEOP  -  Display  of  exterior  orientation  parameter  obs. 

XV_DSP_IMAGES  -  Display  of  image  data 

XV_DSP_C8T  -  Display  of  object  point  cotdinates 

XV_DSP_CAMERA  -  Display  trf  camera  calibraticxi  data 
XV_DSP_SOLEOP  -  Display  of  cc^  results 

XV_DSP_SOLOPC  -  Display  of  results 

XV_DSP_SOLRESID  -  Display  of  residuals 

XV_DSP_SOLNORMAL  -  Display  of  ncamal  equation  system 

6.  Dependencies  (See  Figure  1) 

•  The  module  is  part  of  module  REPORT 

•  The  module  interacts  with  the  modules 

FLOW_CX)N 

CONTR.PAR 

MC»EL_DATA 

8.  Resources  Hardware:  screen  of  a  workstation  or  terminal. 
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8.2  RE_PRINT  -  Output  of  results  and  data  on  printer 

1.  Identification  RE_PRINT 

2.  Type  Module. 

3.  Purpose  Fonnat  and  output  results  and  data  as  reports  and  listings  on  paper  Hottine 

may  be  included  in  an  extended  version  of  die  program. 

Intermediate  results  for  later  resuming  of  the  triangulatkm  is  capable  to  be 
stored  m  files  on  disk. 

4.  Function  Functions  provide  an  interface  to  the  printer  for  output  of  results  and  data. 

5.  Subordinates  This  module  is  ccxnposed  of  the  following  functims: 

XV_PR_FILE 
XV_PR_OIPC 
XV_PR_OOPC 
XV_PR_OEOP 
XV_PR_EC*> 

XV_PR_OPT 
XV_PR_CAMERA 
XV_PR_GENSTAT 
XV_PR_NORMAL 


6.  Dependencies  (See  Figure  1) 

•  The  module  is  part  of  module  REPORT 

•  The  module  intmcts  witii  the  modules 

FLOW_CX)N 

CONTR_PAR 

MOCCLJ5ATA 

8.  Resources  Harware:  laser  printer 

matrix  printer 
(plotter) 


-  I^t  of  a  specific  file 

-  Print  of  ipc  obs.  in  measured  sequence 

-  Prut  rrf  objet^  point  coordinate  obsc^ations 

-  Prmt  of  exterior  (mentation  parameter  obs. 

-  ftmt  of  extaior  orientatitm  parameters 

-  Print  of  object  point  coidinates 
*  Prmt  (tf  camera  calibration  data 

-  Prmt  of  general  statistics 

-  Print  of  normal  equation  system 
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8.3  RE_STORE  -  Storage  of  data  and  results  on  disk 

1.  Identification  RE_STORE 

2.  Type  Module. 

3.  Purpose  Storage  of  input  data  (observations  and  objects),  of  inteimediate  results  and 

the  solutimi  on  disk. 

4.  Function  The  staring  of  die  different  data  will  be  excuted  in  different  functions . 

5 .  Subordiruues  This  module  is  composed  of  die  following  functions: 

ST_CAMERA  -  Saves  camera  calibration  data  in  file 

ST_PARAM  -  Saves  control  parameter  in  file 

ST_RESULT  -  Saves  die  inteimediate  or  final  result  in  file 

ST_EXTORI  -  Saves  exterior  orientation  parameter  in  file 

ST_OBJECT  -Saves  object  point  cooidinates  in  file 

ST_NORMAL  -  Saves  the  ncnnal  equation  system  in  file 

6.  Dependencies  (See  Figure  1) 

•  The  module  is  part  of  module  REPORT 

*  The  module  interacts  widi  the  modules 

FLOW.CON 

E(i.SYST 

CONTR.PAR 

MODEL_DATA 

8.  Resources  Hardware:  disk  few  output 

9.  FUNC_LIB  -  Library  for  common  functions 

1 .  Identification  FUNC_LIB 

2.  Type  Module. 

3.  Purpose  Collection  of  functions  us«i  by  different  modules. 

4.  Function  This  module  contains  functions  for  conqmting  mathemaripai  appiinyrir^is. 

e.g.  rotation  matrix,  design  matrix  coefficients,  and  for  administrating 
common  triplications,  e.g.  common  input  of  values,  file-handling. 

5.  Subordinates  This  module  is  composed  of  the  fdlowing  modules: 

FL.MATH  -  Module  with  mathematical  functions 
FL_UTIL  -  Module  with  utilities 

FL_MATOP  -  Module  with  matrix  operations 

FL_TRAFO  -  Module  with  transformation  functions 
FL.TTME  -  Module  with  time  functions 

6.  Dependencies  (See  Hgure  1) 

•  The  module  is  part  cf  module  (X^TRIS 

*  The  module  interacts  with  the  modules 

INITIAL 

ADJUST 

FLOW.OON 

REPORT 

8.  Resources  '  None  required. 


9.1.  FL_MATH  -  Mathematical  functions 

1 .  Identification  FL3IATH 

2.  Type  Module. 

3.  Purpose  CoUecdon  of  mathematical  functions  used  by  different  modules. 

4.  Function  This  module  omtains  a  fimctkm  for  computing  die  dements  of  die  rotation 

matrix,  the  detenmnant  of  die  rotation  matrix,  an  area  of  an  ot^ue  plane 
triangle  and  an  area  by  the  Gaussian  area  algorithm. 

5.  Subordinates  This  module  is  composed  of  die  following  functions: 

M_ANGRMAT  -  Cooqiutes  the  elements  of  die  rotation  matrix 
M_ANGRMAT_Hn  -  Comp,  the  elem.  of  the  rot  mat.  derived  by  phi 

M_ANGRMAT_(WdE  -  Conqi.  the  elem.  of  the  rot  mat  deriv.  by  ome 

M_ANGRMAT_KAP  -  Comp,  the  elem.  of  the  rot  mat  dmiv.  by  kap 

M_ANGRMAT_TRANS  -  Confutes  transposed  rotation  matrix 

M_ANGRMAT_TRANSPHI  -  Comp,  tranqiosed  rot  mat  deriv.  phi 
M_ANGRMAT_'niANSOME  -  Comp,  transposed  rot  mat  deriv.  ome 
M_ANGRMAT_TRANSKAP  -  Comp,  transposed  rot  mat,  deriv.  kap 
M.ANGLES  -  Ccm^tes  angles  firm  elenomts  of  rot  matrix 
M_MUL1MAT  •  Mult^licates  two  rotation  matrices 

MJ3E1I1MAT  -  Cooqiutes  die  determinant  of  rotation  matrix 
M_AREATRI  -  Conqiiutes  an  area  of  an  oUique  plane  triangle 
M_AREA  -  Coo^utesan  area  by  the  Gaussian  area  algorithm 

6.  Dependencies  •  The  module  is  part  of  module  FUNC_LIB. 

•  Access  from  the  following  modules 

ADJUST 

INITIAL 

8.  Resources  None  required. 
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9.2.  FL_UTIL  -  Utilities  functions 

1 .  Ident^cation  FL_imL 

2.  Type  Module. 

3.  Purpose  Collection  of  utilities  functions  used  by  different  modules. 

4.  Function  This  module  contains  functions  for  input  of  integer,  double  values  and 

words,  also  for  tuning  of  files  for  raiding  and  writing,  checking  die 
presence  of  obse^atims  to  avoid  overwriting  of  data,  finding  observations 
in  the  buffer. 

5.  Subordinates  This  module  is  cmnposed  the  following  functions: 

INP_SKP  -  Cancel  the  rest  cf  an  line 

INP_WRD  -  Input  a  single  word 

INP_INT  -  Input  a  single  integer 

INP.DBL  -  Ii^ut  a  single  double 

FILOPEN.R  -Open  a  file  for  reading 
FILOreN_W  -  Open  a  file  for  writing 
FILOPEN.CAMALL  •  Open  all  files  of  camera  calibration  data 

INCHECKJOEOP  -  Checks  the  presence  of  eop-observatims 

FIND_FOIPC  -  find  the  oipc-observaiion  in  the  buffer 
FIND_FOOPC  -  find  the  o(^-observation  in  the  buffer 
FIND_FOEOP  -  find  the  oe^observation  in  the  buffer 
FIND_IMG  -  find  the  image  in  the  buffer 


<S.  Dependencies 


8.  Resources 


•  The  module  is  part  cf  module 

•  Access  from  fdlowing  modules 

USER.INT 

INPUT 

FLOW.CON 

REPORT 

None  required. 


FUNC_LIB. 


9.3.  FL_MATOP  -  Matrices  operation  functions 

].  Identification  FI._MATOP 

2.  Type  Module. 

3.  Purpose  Collection  of  functions  for  matrices  operation 

4.  Function  Hiis  module  contains  functions  for  printing  and  handling  of  matrices  used 

far3D-Helmerttransfonoation  and  also  for  relative  orientatkm. 

5.  Subordinates  This  module  is  composed  of  the  following  functions: 

ADJ_P11MS  -  Prints  matrix  on  screen 

AD1_UPRMS  -  Prints  iqtper  triangular  matrix  on  screen 

ADIJVfAlMUL  -  Multiplu^  matrices 

ADJ^UATAAHL  -  Com^ites  matrices  sequentially 

ADJ_UADJ2  -  Solves  normals  by  Gtmsaan  factorization 

ADJ_ULIJDEC1  -  LU  deconqxrsition  of  upper  triangular  matrix 

ADJ_ULUSOLl  -  Sieves  LU  decomposed  system 

ADJ_UMA'i'ZltkO  -  Sets  all  elements  m  iqrper  triang.  matrix  to  zero 

ADJ_MATZER0  -.Sets  all  elements  of  a  ^  matrix  to  zero 

ROTMAT  -  Computes  rotation  matrix  (different  way) 

FC81M..ABC  -  Forms  partials 

VEC_PROD  -  Mult^licates  vectors 

VECJFCXIM  -  Changes  storage  of  matrices 

MAT_INV_PIV  -  Invertes  a  sym.  fully  stared  matrix  (pivoting) 

FACT_BACK_INV  -  Factorizes,  substitutes  and  invertes  normals 

6.  Dependencies  •  The  module  is  part  of  module  FUNC_LIB. 

*  Access  Burn  fallowing  modules 

INITIAL 

FUNC_LIB 

8.  Resources  None  required. 

9.4.  FL_TRAFO  •  Library  of  functions  for  transformation 

1.  Identification  FL_TRAFO 

2.  Type  Module. 

3.  Purpose  Collection  of  functions  for  SD-Helmert  transformation. 

4.  Function  This  module  contains  functions  for  pttmaring  the  data  for  the  transformation, 

for  computing  the  partial  derivatives  tx  each  observation  and  ftx-  computing 
the  seven  parameter  of  transformation  an  adjustment 

5.  Subordinates  This  module  is  composed  of  the  fcrilowing  functions: 

T_AI3PAR  -  I^ta  preparation  for  transfonnation 

TJIRAPO  -  3D-Helmen  transformation 

T_PAIt_OPC  -  Partials  of  each  observation 

6.  Dependencies  (See  Figure  1) 

*  The  module  is  pan  of  module  FUNC.UB. 

» Access  firom  the  modules 

FLOW.OON 

8.  Resources  None  required. 


9.5.  FL_TIME  -  Time  fkmctions 

1.  Ident^cation  FL.TIME 

2.  Typ€  Module. 

3.  Purpose  Collection  of  functions  for  measuring  the  cpu,ekq>sed  time. 

4.  Function  This  module  contains  fiinctions  to  measure  the  cpu  andlAir  eitq>sed  time,  to 

return  the  cunrat  date. 

5.  Subordinates  This  module  is  conqrosed  of  the  following  functions: 

TIMER^START_CPU  -  Start  tire  dock 

TIMEILSTOP.CPU  •  Stop  the  clock  and  measure  cpu-time 
HMEILSTART  -  Start  tiie  dock 

TIMEIL.STOP  -  Stop  the  dock  and  measure  the  dt^sed  time 

GET_CURRENT_DATE  JNT  -  Get  the  current  date  as  integer 

GET_CURRENT_DATE_CHR  -  Get  the  current  date  as  char 

6.  Dependencies  •  The  module  is  part  of  module  FUNC.LEB. 

•  Access  firm  the  following  modules 

ADJUST 

8.  Resources  None  required. 

10.  MEM_ALLOC  •  Memory  allocation  functions 

1.  Identification  hJOSM^JiUClC 

2.  Type  Module. 

3.  Purpose  Cdlection  of  fiuictions  for  memosy  allocation  and  teallocation  for  specific 

items  of  the  data  structures  defined  in  the  system. 

4.  Function  This  module  contains  functions  which  aDocaie  and  reallocate  space  for 

objects  of  specified  types,  defined  in  this  system. 

5.  Subordinates  This  module  is  con^osed  of  die  following  modules: 

MEM_<%S  -hfemoiy  allocation  for  observation  smictuies 

MEMJ*AR  -Menxxy  allocation  for  parameter  structures 

MEMJSQSYS  -Memory  allocation  for  normal  equation  structures 

6.  Dependencies  (See  Figure  1) 

•  The  module  is  part  of  modide  OLTRIS 

•  The  module  interacts  with  die  noodules 

MODEU3ATA 

EQ_SYST 

8.  Resources  liatdwae:  -  REM  (Random  excess  memory). 


10.1.  MEM_OBS  •  Memory  allocation  for  observation  structures 

1 .  Identyication 

2.  Type  Module. 

3.  Purpose  Ccrilection  of  functioas  for  memory  allocatkn  and  release  for  dw  data 

structures  of  observations,  defined  in  the  system. 

4.  Function  Ttiis  module  contains  functions  which  allocate  and  release  ^>ace  for 

observations  of  qtecified  types,  defined  in  this  system. 

5.  SubortUncaes  This  module  is  cott^tosed  of  the  following  functions: 

MEM_OIPC  -Monoty  allocatioo  for  tm^pomt  observations 

MEM-OQPC  -Memory  allocation  for  objm  point  observations 

MEM_C£C^  -  Menxsy  allocation  for  ext  orient  parano.  observ. 

MEMF_(HPC  -Releaseofmenoory  for  im^  point  obsetvarx»s 
MEMF_OCK*C  -Rekaseof  memory  for  objtM  point  observations 
MEMFjC^CX*  -  Release  of  memory  for  ext  orient  parattt  observ. 

6.  Dependencies  (See  Figure  1) 

•  The  module  is  part  of  module  MEM_ALLOC 

•  The  module  interacts  with  the  modules 

MODEL_PATA 

INPUT 

8.  Resources  Hardware:  -  REM  (Random  excess  memory). 
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10.2.  MEM_PAR  •  Memory  allocation  for  parameter  structures 

1 .  Identification  MEM_PAR 

2.  Type  Moduk. 

3.  Purpose  CMlectkMoffii«aioosformemoiy  allocation  and  release  of  metnoiy  for 

data  stiuctuies  of  parameters  defined  in  the  systeia 

4.  Function  This  module  contains  functions  which  aOocate  and  release  space  for 

parameters  of  q)ecified  types,  defined  in  this  system. 

5.  Subordinates  This  module  is  coiiq>osed  of  the  following  functions: 

MEM_OPT  -Memory  allocation  for  object  point  item 

MEM_CX’C  -Memory  allocation  for  otgect  point  coordinates 

M^_EOP  -  Memory  allocation  forest  orient 
MEM_IMG  -  Memory  allocatioo  for  observation  structures 

MEM_CAMERA  -  Memory  allocation  for  parameter  structures 

MEMF_OPT  -  Release  of  memory  for  object  point  item 

MEMF_OPC  -  Release  of  memory  for  otject  point  coordinates 
MEMFJECX’  -  Relrase  of  memory  for  ext  orient 

MEMF_IMG  -  Release  of  memory  for  observation  structures 

MEMF_CAMERA  -  Release  of  memory  fm^  parameter  structures 

6.  Dependencies  (See  Figure  1) 

•  The  module  is  part  of  module  MEM.ALLjOC 

•  The  module  interacts  widi  the  modules 

MODEL_DATA 

INPUT 

8.  Resources  Hardware:  -  REM  (Random  excess  memory). 


10.3.  MEM_EQSYS  •  Mem.  alloc,  for  normal  equation  structures 

1.  Identification  MEM^QSYS 

2.  Type  Module. 

3.  Purpose  Ccdkction  of  function  for  meinoiyallocatioa,reallocatun  and  idease  of 

memcry  for  q)ecific  items  of  the  normal  equatum  data  structures  defined  in 
the  system. 

4.  Function  This  module  oontaims  fanetjons  which  allncatf.,  rpaUnratf.  fliy)  lylewsr  spare 

for  objects  of  normal  equation  system,  defined  in  diis  system. 

5.  Subordinates  This  module  is  composed  of  the  following  functions: 

MEM_UFPER  -  Memory  aUocatioa  for  upper  triangle  of  matrix 

MEM_VECl<jR  -hfonoory  allocation  for  c^umn  or  row  vector 
MEMLAIPC  -Memory  allocation  for  design  submatiix  of  q>c 
MEM..ACX*C  -  Memory  allocation  fm  design  submatiix  cf  opc 

MEM^AEOP  -Memory  allocation  for  design  submatiix  of  eop 

MEM_UXT  -  Memory  aUocatioa  for  submatrix  in  the  hypenow 
MEM_HYP  -  Memory  allocation  for  hyperline  in  the  nmmals 

MEMR_UPreR  -Memory  reallocation  for  upper  tiian^rf  matrix 
MEMR_VECrCXt  -  Memory  reallocation  for  c^umn  or  row  vector 

MEMF.UFPER  -  Release  of  memory  for  upper  triangle  of  matrix 
MEMF_VECrOR  -Releaseofmemory  for  column  or  row  vector 
MEMF.^AIPC  -  Release  of  memory  for  design  submatiix  of  ipc 

MEMF_AOTC  -  Release  of  memory  for  design  submatiix  of  opc 

MEMF_AE(^  -  Release  of  memory  for  design  subaiatrix  of  cop 

MEMF_UXT  -  Release  rtf  memory  for  submatiix  in  tiic  hypenow 

MEMF_HYP  -Releaseofmemory  for  hyperline  in  the  nramals 

6.  Dependencies  (See  Hgure  1) 

•  The  module  is  part  of  module  MEM-ALLOC 

•  The  module  interacts  widi  the  modules 

MCXJEL^DATA 

EQ_SYST 

8.  Resources  Hardware:  -  REM  (Random  excess  memory). 
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11.  CONTR  PAR  -  Control  parameters 

1.  Identification  CX)NTRJ*AR 

2.  Type  Data  module. 

3.  Purpose  To  contain  centred  panuneteis,i.e.  all  parameteis  and  data  used  for 

controlling  die  different  fiinctkns  performed  by  the  system. 

4.  Function  The  module  contains  the  initialization  of  control  flow  parameters,  as  well  as 

critical  values  and  other  general  parameters  needed  by  some  computations. 

5.  Subordinates  This  modute  contains  the  following  functions; 

CP_ADJUST  -  Control  parameters  for  the  adjustment 
CP_EQSYS  -  Cratnd  parameters  for  the  normal  equation  system 
CP.hdSC  -  Control  parameters  for  the  miscellaneous  plications 

6.  Dependencies  *  This  module  is  part  of  module  OLTRIS 

•  Interacts  probaldy  with  all  high-level  functional  modules. 

8.  Resources  The  data  is  supposed  to  reside  in  core. 

12.  MODEL_DATA  >  Data  for  entities  of  the  mathematical 
model 

1.  Identification  M(X)EL_DATA 

2.  Type  Data  module. 

3.  Purpose  Orgamre  storage  of  data  concerning  entities  used  in  die  mathematical  model. 

4.  Function  This  module  contains  data  structures  for  the  "physical"  entities  cameras, 

inages,  and  objea  p^ts;  for  a  priori  and  conqwted  parameters,  and  for 
obsmatkms  and  design  matrix. 

5.  Subordinates  This  module  is  composed  of  the  modules  (see  figure  4): 


C»>T 

-Object  point  data 

OPC 

-  Ot^ect  pc  lUt  cootdinaies 

EC»> 

-  Exterior  ciientation  parameters 

IMAGE 

-  Image 

CAMERA 

-Oumadata 

OIPC 

-  Observations  of  image  point  coordinates 

OOPC 

-  Observations  of  object  point  coordinates 

OEOP 

-  Observations  of  exterior  orientation  parameters 

MD.FUNC 

-  Functions  for  handling  the  data  structures 

6.  Dqjendencies  •  This  mndnle  is  pan  of  module  (X^HUS 

•  Supposed  to  be  accessed  from  the  modules 

INPUT 

INITIAL 

ADJUST 

QUAL_CON 

REP(»T 

•  Uses 

MEhfjVLLOC 

•  By  pointeis  and  indices  cioss-ielated  to  data  in 

EQ_SYST 

8.  Resources  The  data  is  supposed  to  reside  in  core. 


Figure  4:  Reladons  between  ifae  dna  modnles  of  MOKL_DATA  Ihd  or  me  Htowbeads  indicate  whedier 
each  m«tnnce  nf  die  entity  M  lha  nther  wid  nf  ihe  Hne  may  he  irtowl  m  leveial  nr  only  one  inaance  of  the 
entity  pointed  at.  The  greyed  boxes  belong  to  module  EQ.SYST. 
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12.1.  MD_FUNC  -  Functions  for  handling  the  data  structures 

1 .  Ideni^cation  MD_FUNC 

2.  Type  Module. 

3.  Purpose  Organize  the  handling  of  the  ooodeLdata  stroctures. 

4.  Function  This  module  contains  functions  for  ii^ut  ct  data  into  the  modd.data 

structures  and  for  finding  data  in  its  lists. 

5.  Subordinates  This  module  is  conqx>sed  of  the  following  functions: 

MDLOlPC  -  Iiq>ut  obsetvatkns  of  im^e  point  coordinates  into  list 

MDLOOPC  -  Input  observations  of  objM  point  coordinates  into  list 

MDI_OEOP  -  Input  observations  of  exL  oiientation  param.  into  list 
MDLCX’T  -  Input  object  point  coordinates  into  list 

MDL  A(G  -  Ii4>ut  image  data  into  list 

MDI_CAMERA  -  Ii^ut  camera  data  into  list 

MDF_OIPC  -  Find  observations  of  image  point  coordinates  in  list 

MDF.OCX’C  -  Find  observations  of  object  point  cotHdinates  in  list 

MDF_OG(^  -  I^nd  observadcms  of  ext  orientation  parameters  in  list 

MDF.CXn*  -Find  object  point  comdinates  in  list 

MDF_!^G  -  Find  image  data  in  list 
MDFjCAMERA  -  Hnd  camera  data  in  list 

6.  Dependencies  *  This  module  is  part  of  module  MODEL.DATA. 

*  Access  fiom  die  modules 

INPUT 

INITIAL 

ADJUST 

QUAL_CON 

REPORT 

•  Uses 

data  modules  of  MO[M._DATA 
8.  Resources  None  required. 


13.  EQ_SYST  -  Normal  equation  system  and  related  data 

1 .  Identification  EQ_SYST 

2.  Type  Data  module. 

3.  Purpose  Organize  storage  for  the  normal  equation  system. 

4.  Function  This  module  cmttains  data  conqnising  die  nonnal  equatkm  system,  the 

elements  of  which  are  co^Bdent  matrix,  constant  vectotfs)  and  soludm 
vector.  Also  included  is  related  woik-space,  e.g.  for  obs^admi  equations 
during  sequential  iqidatts. 

5.  Subordinates  This  module  is  conqiosed  (rf  the  modules: 

EQJB  ASB  -  Basic  parts  of  the  normal  equations  ^stem 

EQ_HYI^  •  Hyper-row  parts  of  the  nor^  equations  system 
EQ_WORKSP  -Woic-space  related  to  the  nonnalequaticxis  system 

EQ_FUNC  -  Functions  for  handling  die  nmnal  equations  system 

6.  Dependencies  •  This  module  is  part  of  module  OLTRIS 

•  Access  fiom  the  modules 

ADJUST 

QUAL_CON 

•  Uses 

MENLALLOC 

8.  Resources  To  achieve  minimal  re^xmse  time,  this  data  should  reside  permaiiendy  in 
core.  If  possiUe  it  shoidd  be  prevented  Bom  any  swapping  /  paging. 

13.1.  EQ  BASE  •  Basic  parts  of  the  normal  equations  system 

1 .  Identification  EQBASE 

2.  Type  Data  module. 

3.  Purpose  Organize  storage  for  the  basic  parts  of  die  normal  equation  system. 

4.  Function  This  module  contains  data  cooqxising  the  basic  pan  of  the  normal  equation 

system,  the  elements  of  sdiich  are  coefiicient  nuurix,  constant  vectorfs)  and 
soluticHi  vector. 

5.  Subordifuites  This  module  is  conqiosed  of  the  ftdlowing  data  structures: 

A_IPC  -  Design  submatrix  for  im^  point  observation 

A_OPC  -  Design  submatrix  for  ob)^  point  observation 

A_EOP  -  Design  submatrix  for  ext  orient  paratn.  observation 

6.  Dependeraies  •  This  module  is  pan  of  module  EQ_S  YST. 

•  Access  Bom  the  modules 

ADJUST 

(yjAL_CXMf 

•Uses 

MEMjMJjOC 

8.  Resources  To  achieve  minimal  remonse  time,  this  data  should  reside  permat^dy  in 
core.  If  possible  it  should  be  prevented  Bom  any  swqifnng  /  paging. 
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13.2.  EQ_HYPER  -  Hyperrow  parts  of  the  normal  equat.  system 

1 .  Identification  E{i_HYPER 

2.  Type  Data  module. 

3.  Purpose  Organize  storage  for  the  hypeirow  parts  of  the  normal  equation  system. 

4.  Function  This  module  contains  data  comprising  the  hyperrow  parts  of  the  nonnal 

equaticm  system,  the  elements  of  which  are  coefficient  matrix,  constant 
vector(s)  and  solution  vector. 

5.  Subordinates  This  module  is  composed  of  the  data  structures: 

HYP_OPC  -  Hypcriine  in  the  normals 

SUB_XT  -  Submatrix  in  the  object  point  hyperrow 

6.  Dependencies  •  This  module  is  part  of  module  EQ_S  YST. 

•  Access  from  the  modules 

ADJUST 

QUAL_CON 

•  Uses 

MEM_ALLOC 

8.  Resources  To  achieve  minimal  response  time,  this  data  should  reside  permanently  in 
core.  If  possible  it  shoidd  be  prevented  from  any  swapping  /  paging. 

13.3.  EQ_WORKSP  •  Workspace  related  to  the  normal  eq.  syst. 

1 .  Identification  EQ_WORKSP 

2.  Type  Data  module. 

3.  Purpose  Organize  storage  of  related  dataA'ariables  for  the  normal  equation  system. 

4.  Function  This  module  contains  data  coraimsing  related  dataA'ariables  of  the  ncnmal 

equation  system.  Also  included  is  allocated  woik-space  of  this  gobal 
variables. 

5.  Subordinates  This  module  is  composed  of  die  following  data  (variables): 

IMAGE_INIT  -  Initia]  number  of  images  for  allocation 
IMAGE_AIX)  -  Number  of  additional  images  at  reallocation 
E(i_PARAMS  -Number  of  image  parameter  per  image 
N_AWSP  -  Number  of  lines  (A)s.  eq.  workspace  to  allocate 

SEZE_HYP  -  Number  of  lines  to  allocate  in  a  hyperrow 
EQ_IND  -  Compute  linear'  index  in  an  'Upper'  matrix 

6.  Dependencies  •  This  module  is  part  of  module  E(i.SYST. 

•  Access  from  the  modules 

ADJUST 

QUAL_CON 

8.  Resources  To  achieve  minimal  response  time,  this  data  should  reside  permanently  in 

cote.  If  possible  it  should  be  prevented  from  any  swapfug  /  paging. 
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13.4.  EQ_FUNC  •  Functions  for  handling  the  normal  eq.  system 

1.  Identification  EQ_FUNC 

2.  Type  Module. 

3.  Purpose  Orsanize  administrating  of  the  normal  equation  system. 

4.  Function  This  module  contains  functions  fiv  administTating  and  house-keq>ing  of  the 

n<»mal  equation  system. 

5.  Subordinates  This  module  is  composed  of  tiie  following  functions: 

EQt_INIT  -  Initialize  the  normal  equation  syston  storage 

EQ_EXPAND  -  Expand  the  normals 
EQ_CLEAR  -  Clear  the  normal  equations  sys^n 

EQ.SETHYP  -  Basic  parts  of  the  normal  equations  system 

EQ_2ERO  -  Zeio-ffl  tiie  cunent  normal  equations  system 

EQ_ntINT  -  Print  the  current  nonnal  equations  system 
EQ_S17FAL  -  Return  sizes  o(  currently  allocatcxl  q>ace  fen-  normals 
EQ_SIZCU  -  Return  sizes  of  the  current  normal  equations  system 
EQ_POINTR  -  Return  pmnters  to  ntmnal  equation  system  storage 

EQ_HYITjOAO  -  Load  a  hyperrow  frmn  separate  storage 
EQ_HYPS  AVE  -  Save  a  hyperrow  to  separate  storage 

6.  Dependencies  •  This  module  is  part  of  module  EQ_SYST. 

•  Interacts  with  the  following  modules 

ADJUST 

QUAL.CON 

S.  Resources  None  requited. 


